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ABSTRACT 

By  a  re-definition  of  the  standard  coefficients  of  form, 
wherein  displacement  is  weighted  by  its  distribution  in  depth, 
the  analysis  of  wave  making  resistance  is  simplified.   Using 
the  Taylor  Standard  Series  as  a  data  base,  regression  analyses 
are  employed  to  develop  a  universal  wave  drag  function  of  only 
one  variable  for  each  speed  length  ratio. 

Form  drag  is  decomposed  into  components  which  facilitate 
analysis.  Prediction  formulae  are  developed  graphically  from 
the  Taylor  Standard  Series  data. 

Using  the  prediction  formulae  developed,  wave  drag  and 
form  drag  are  predicted  for  a  number  of  sample  ships  and  are 
compared  to  that  predicted  by  Gertler's  reanalysis  of  the 
Taylor  Standard  Series.   Similar  comparisons  are  made  for  a 
half  dozen  models  which  are  unlike  the  Taylor  Standard  Series 
parent  with  their  individual  towing  tank  test  results.   Dif- 
ferences and  possible  sources  of  error  are  discussed. 
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NOMENCLATURE 

A  sectional  area 

A*  sectional  area  weighted  by  its  distribution  in  depth 

A  midship  section  area 

A*  midship  section  area  weighted  by  its  distribution  in 

m  depth 

A[0]  free  wave  amplitude  spectrum 

B  beam  at  midship  section 

B/T  beam  to  draft  ratio 

C_,  block  coefficient  =  V/LBT 

Cf  form  drag  coefficient  =  R-/1/2  pV2S 

C  skin  friction  error  coefficient  =  R  /1/2  pV2L2 
o  o 

C,  coefficient  of  the  form  drag  coefficient  proportional 


'1 


2T  2 


to  displacement  =  R,/l/2  pV  L 


C  midship  section  coefficient  =  A  /BT 
m  .  *  nr 

C*  midship  section  coefficient  of  the  weighted  sectional 

m  area  ^  A*/BT 

nr 

C  longitudinal  prismatic  coefficient  =  V/A  L 

C*  weighted  longitudinal  prismatic  coefficient  =  V*/A*L 

C  midshio  section  coefficient  =  A  /BT 
x  x 

CR  residuary  resistance  coefficient  =  R  /1/2  pV2S 

C_(.5)  residuary  resistance  coefficient  at  V//L  =  .5 

C„  wetted  surface  coefficient  =  S//VL 


C       volumetric  coefficient  =  V/L 
v 


C*      weighted  volumetric  coefficient  =  V*/L 


2T  2 


C       wave  drag  coefficient  =  R  /1/2  pV  L 

w  3  w 
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V  displaced  volume 

V*  displaced  volume  weighted  by  its  distribution  in  depth 

5  boundary  layer  thickness 

F  Froude  number  =  V//gL,  where  V  is  in  feet  per  second 

F  skin  friction  error  function 
o 

F,  form  drag  slope  function 

2 

Fp  interim  wave  drag  function  =  (F.J  (1000) * 

F..  universal  wave  drag  function 

g  gravitational  constant  =  32.2  ft/sec2 

K  wave  number  =  2tt/A  =  g/V2 

L  length  on  the  waterline 

LCB  longitudinal  center  of  buoyancy 

LCB*  longitudinal  center  of  the  weighted  displaced  volume, 

V* 

L/T  length  to  draft  ratio 

A  wave  length 

<J)  velocity  potential  of  a  2-D  wave 

cf>  velocity  potential  of  a  2-D  wave  at  z  =  0 

Rf  form  drag 

R  residuary  resistance 

R  wave  drag 

w  r 

R  skin  friction  resistance  error 
o 

R,  form  drag  component  proportional  to  displacement 

S  wetted  surface 

T  draft  at  the  midship  section 
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0  direction  of  propagation  relative  to  ship  heading 

V//L  speed  length  ratio  (L  in  ft.,  V  in  kts.) 

V  ship  speed  (in  knots  or  ft.  per  second) 

y(x,z)  offset  of  hull  form 

y  half  width  of  separated  wake 
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I .   INTRODUCTION 

Although  many  phenomena  encountered  in  Naval  Architecture 
have  been  reduced  to  relatively  simple  prediction  formulae 
with  reasonable  engineering  accuracy,  there  are  as  many  which 
do  not  admit  to  any  simple  analysis.   Among  the  most  notorious 
are  those  phenomena  which  are  responsible  for  the  residuary 
resistance  of  ships.   Withstanding  centuries  of  scrutiny  and 
data  taking,  they  have  steadfastly  refused  to  collapse  into 
any  reasonably  simple  relationships  involving  the  standard 
coefficients  of  ship  hull  form. 

During  the  18th  century,  Chapman  [3]  built  ship  model 
towing  tanks  and  developed  scaling  laws  and  ship  synthesis 
models  from  his  observations.   Unable  to  sort  out  all  the 
phenomena  involved  in  ship  resistance  with  his  ship  model  ex- 
periments and  other  observations,  he  complained,  "For  ships, 
we  have  to  fear  an  infinity  of  bad  qualities  of  the  greatest 
consequence,  which  we  are  never  sure  of  being  able  to  remove, 
without  understanding  the  theory." 

During  the  last  century  since  the  publication  of  Froude ' s 
paper  [13]  in  1868,  ship  model  basins  have  been  constructed  in 
many  parts  of  the  world  for  the  purpose  of  measuring  the  re- 
siduary component  of  ship  resistance. 

Although  we  have  probably  sorted  out  most  of  the  phenomena 
involved  in  ship  resistance  and  have  explained  them  with  the- 
ories, we  have  not  yet  succeeded  in  predicting  residuary  drag 
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in  terms  of  standard  ship  form  coefficients.   This  is  due  not 
only  to  the  mathematical  complexity  of  the  theory,  but  also  to 
the  quantities  required  to  be  measured  by  the  theory  which  are 
in  general  not  those  measured  by  standard  coefficients  of  ship 
hull  form.   The  difficulty  is  best  expressed  by  Weinblum  [26] . 
"...  the  wave  resistance  depends,  to  a  first  approximation, 
upon  a  complicated  function  of  the  surface  slope  in  the  longi- 
tudinal direction,  i.e.,  on  derivatives.   On  the  other  hand, 
the  most  commonly  used  hull  coefficients  are  integrals  which, 
even  when  kept  constant,  still  admit  of  very  wide  variations 
of  the  slopes.   We  realize  now  why  the  solution  of  the  basic 
problem  of  the  model  basins  mentioned  above  -  to  establish  the 
resistance  as  a  function  of  the  form  -  remains  almost  hopeless 
as  long  as  the  ship  surfaces  (or  at  least  their  most  important 
features)  are  not  defined  in  a  rigorous  way  by  mathematical 
expressions. " 

It  is  the  purpose  of  this  thesis  to  do  just  that,  i.e., 
to  relate  residuary  drag  to  hull  form.   The  success  of  this 
endeavor  depends  on  two  considerations.   The  first  is  that  the 
description  of  the  hull  form  adopted  for  the  thesis  contains 
parameters  which  measure  the  distribution  of  displacement  in 
depth  and  its  effect  on  the  velocity  potential  with  Froude 
number,  in  place  of  some  of  the  standard  coefficients  of  form. 
The  second  is  the  word  "almost"  in  Weinblum' s  quote. 
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II.   WAVE  DRAG 

The  wave  drag  expressed  in  terms  of  the  generated  free 
wave  amplitude  spectrum,  A(9),  will  be: 


1      2 
R  =  i  TipV2 
w    2   ^ 


+  7T/2 

f 

[A(6) ] 2cos36de  (1) 

-tt/2 


The  cos30  terni  represents  the  weighting  of  the  contribution  of 
each  part  of  the  wave  spectrum  to  wave  drag  with  the  relative 
direction  of  propagation.   The  cos30  weighting  function  indi- 
cates that  the  major  contribution  to  wave  drag  will  be  from 
those  waves  travelling  in  the  same  direction  as  the  ship  (9  =  0). 
A  Kelvin  wave  pattern  is  illustrated  in  Figure  1  (taken  from 
reference  22)  which  shows  ship  generated  waves  propagating  in 
various  directions.   In  this  illustration,  it  is  the  transverse 
waves  following  the  ship  which  radiate  most  of  the  energy  away 
from  the  ship.   In  view  of  the  relatively  small  contribution 
to  wave  drag  from  the  shorter  diverging  waves,  it  is  a  reason- 
able simplification  to  represent  the  wave  system  causing  wave 
drag  by  the  wave  length  and  wave  number  of  the  predominate 
transverse  wave. 

Several  analytical  methods  have  been  derived  predicting 
the  wave  drag  of  slender  bodies.   Mitchell's  [14]  integral  is 
probably  the  method  most  often  cited.   Rather  than  standard 
coefficients  of  form,  Mitchell's  integral  measures  the  form  of 
the  hull  in  terms  of  longitudinal  slopes  of  the  surface  of  the 


<v.  » 


-First  14  crests  of  a  Kelvin  wave  group  caused 
by  a  travelling  disturbance  at  O. 
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hull  at  each  point  on  the  hull.   The  pressure  change  due  to 
the  resulting  wave  system  is  integrated  over  the  ship  to  de- 
termine the  total  drag  force.   The  coordinate  system  is  shown 
in  Figure  2. 

Using  Mitchell's  integral,  wave  drag  may  be  expressed  as 
(taken  from  reference  26): 


R  =  -2 


<Sp  dy  dz  =  -2 


dp  (-^-)dxdz 


=  -2pV 


Jj 


<H>  <&>*<* 


(2) 


or 


=  tea*. 

ttV* 


(T2  +  J2) 


/X^T 


dX 


(3) 


where 


and 


+e  K(x) 


I  = 


(^)e 
^x; 


A2gz 
V' 


cos  (-r^-)dxdz 


-e  0 


J  = 


+e  K(x)       A2gz 

(|J)a"  V2 

8x 

-e  0 


sin(-r^)dxdz 


where  z  =  K(x)  is  the  profile  of  the  fairbody. 

It  may  be  seen  from  the  foregoing  that  Mitchell's  integral 
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K 

I 

v5 


CX 
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integrates  the  slope  of  a  given  waterline,  -^,  weighted  by 
factor, 

A2qz       A2 


e       -  e  (4) 


where  F  =  — —  is  a  local  Froude  number  based  on  the  ship  speed 
and  the  depth  of  the  waterline,  over  the  length  of  the  ship. 
The  area  of  each  waterline  is  then  weighted  by  its  speed  and 
depth  below  the  free  surface.   The  weighted  waterline  areas 
are  in  turn  integrated  over  the  draft,  which  weights  each  in- 
crement of  displacement  by  its  depth. 

The  "effective"  displacement  of  a  ship  which  contributes 
to  the  wave  making  velocity  potential  is  therefore  something 
less  than  the  real  displacement  by  the  weighting  factor  in 
Mitchell's  integral.   If  the  lines  of  two  vessels  of  equal 
real  displacement  are  similar  except  for  proportional  changes 
in  beam  to  draft  ratio,  the  smaller  beam  to  draft  ratio  should 
have  less  effective  displacement  and  less  wave  drag  than  the 
vessel  with  the  larger  beam  to  draft  ratio.   The  difference  in 
the  effect  of  the  displacement  at  various  depths  is  illustrated 
in  Figure  3 . 

In  addition,  due  to  the  wave  number,  K,  in  the  exponent, 
it  is  expected  that  the  magnitude  of  any  effect  on  wave  drag 
due  to  the  distribution  of  displacement  in  depth  will  also 
depend  on  speed  length  ratio  or  Froude  number.   For  high  speeds 
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where  the  wave  length  of  the  ship  generated  waves  is  approach- 
ing the  waterline  length  of  the  ship,  the  change  in  velocity 
potential  over  the  draft  of  the  ship  will  become  relatively 
small  for  ships  with  large  length  to  draft  ratios. 

The  variation  of  the  velocity  potential  with  depth  may  be 
written  as: 


<|>(z)  =  <J>(z=0)e  KZ  =  <f>oe~KZ  (5) 


where  the  wave  number,  K  =  ^-  =  -y-   and  A  is  the  wavelength  of 
the  transverse  waves  generated  by  the  ship. 
At  high  speed  length  ratios , 
X  approaches  L,  and 
K  approaches  — . 


The  speed  length  ratio, 

V (knots)         .  „  V (knots) 
— - — -== — -  approaches  -~= 


1.689/27rVVg  ~  1'34 


and  the  Froude  number,  F 

V  V 

approaches 


7gT  "ft—"—  /gA 

V  1 


/ZtT 


=  .399 


In  addition,  the  effect  on  the  velocity  potential, 
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2ttz 

— Kz  T 

(J)  e   "  approaches  $   e 


At  this  high  speed,  the  effect  on  the  velocity  potential 
at  various  values  of  z  is: 


<J)(z=j)    =   <f>oe    7T/2   =    .208<|>  (6) 


.  ,  L .  .       -2tt/9.065         1, 

<J,(Z=9T065)    =    V  =    2*0  (7) 


d>(z*^)  =  4>oe  27T/5°  =  .8824>o  (8) 


The  speed  of  a  displacement  vessel  corresponding  to  a 
Froude  number  of  -rw—  is  often  referred  to  as  the  "hull  speed" 
for  that  vessel.   It  may  be  seen  from  the  foregoing  calculations 
that  the  magnitude  of  the  velocity  potential  will  decay  in 
depth  from  the  free  surface  at  a  rate  which  will  reduce  its 
value  to  one  half  that  at  the  surface  at  a  depth  of  about  one 
ninth  of  the  waterline  length  of  the  vessel  at  hull  speed. 

As  ship  speed  is  decreased  from  a  speed  length  ratio  of 
1.34,  i.e.,  down  from  "hull  speed",  wave  drag  decreases  sharp- 
ly.  The  wave  number,  however,  increases  rapidly  as  speed  is 
decreased,  thereby  increasing  the  rate  of  decay  of  the  expo- 
nential with  depth.   There  exists,  therefore,  a  speed  range 
wherein  the  wave  drag  has  not  yet  become  negligibly  small,  but 
where  the  rate  of  decay  with  depth  of  the  velocity  potential 
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has  become  very  strong.   In  this  situation  the  effect  of  the 
distribution  of  displacement  in  depth  on  the  wave  drag  may 
become  very  important.   Figure  4  illustrates  the  change  in  the 
decay  profile  of  the  velocity  potential  with  ship  speed. 

The  decay  of  the  velocity  potential  with  wave  number  and 
depth  may  be  written: 


4.  /  ^       ~Kz    ,   -gz/V2  ._. 

<J)(z)  =  <f>Qe    =  <f>oe  y  '  (9) 


To  demonstrate  the  variation  in  the  decay  curve  of  the 
effect  of  displacement  with  depth,  caused  by  length  draft 
ratio,  we  may  look  at  the  amount  by  which  the  effective  dis- 
placement is  reduced  at  the  keel.   Taking  z=T, 


+  (»  =  *  e"9T/v2   =  *  e-(gVV2)(T/L) 


.  t    e-l/(F2)(L/T) 


where  F  is  the  Froude  number  and  L/T  is  the  length  to  draft 
ratio. 

It  may  be  seen  from  this  approximation  that  a  large  length 
to  draft  ratio  will  flatten  out  the  decay  curve  in  a  manner 
similar  to  that  of  a  large  Froude  number.   Long  shallow  ships 
or  high  speeds  will  therefore  result  in  the  least  influence  of 
the  distribution  of  displacement  over  depth  on  wave  drag.   For 
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either  long  shallow  ships  or  high  speeds,  in  other  words,  the 
effective  displacement  will  approach  the  real  displacement  in 
its  contribution  to  wave  making. 
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III.   FORM  DRAG 

Without  viscosity,  the  water  is  allowed  to  slip  past  the 
hull  of  a  ship,  its  interaction  with  the  ship  being  limited  to 
normal  forces  and  wave  radiation.   D'Alembert's  paradox  allows 
the  pressure  distribution  aft  to  produce  the  same  normal  forces 
as  those  acting  against  the  forward  part  of  the  ship.   With 
the  introduction  of  even  the  smallest  amount  of  viscosity, 
however,  a  no-slip  condition  is  imposed  upon  the  skin  of  the 
ship.   This  would  directly  result  in  the  shearing  of  viscous 
layers  of  fluid  very  close  to  the  hull  and  the  diffusion  of 
vorticity   outward   at   a   rate   proportional   to   the 
gradient   of   velocity   at   the   hull.    In   addition,   at 
high   Reynold's   numbers,   vorticity   is   physically   con- 
vected  away  from  the  hull  by  turbulent  mixing.    By  these 
mechanisms,  momentum  is  exchanged  with  the  flow  and  the  velo- 
city distribution  is  changed  from  that  of  the  inviscid  flow. 
For  very  high  Reynold's  numbers  this  region  of  high  vorticity 
flux  and  momentum  exchange  is  limited  to  a  relatively  thin 
boundary  layer  which  grows  from  the  bow  to  the  stern  and  con- 
tinues aft  as  all  or  part  of  the  wake.   Figure  5  shows  the 
boundary  layer  thickness  and  velocity  profiles  along  the  hull 
of  a  ship.   The  presence  of  the  flow  entrained  in  the  boundary 
layer  prevents  the  streamlines  aft  from  closing  as  quickly  as 
those  in  the  inviscid  case.   As  a  result,  the  pressure  is  not 
able  to  increase  aft  back  up  to  the  value  which  would  produce 
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normal  forces  aft  with  forward  components  to  equal  the  after 
components  of  the  normal  forces  on  the  forward  part  of  the 
ship.   The  difference  in  longitudinal  components  of  the  normal 
forces  resulting  from  a  reduction  in  the  pressure  field  aft  by 
the  boundary  layer  is  the  viscous  pressure  drag. 

If  the  curvature  of  the  hull  lines  aft  is  too  abrupt  or 
if  for  any  other  reason  the  pressure  gradient  in  the  direction 
of  flow  is  strongly  adverse,  the  boundary  layer  flow  may  sepa- 
rate from  the  hull,  leaving  a  region  of  reversed  flow  following 
the  ship.  The  region  of  flow  following  the  ship  constitutes  a 
momentum  defect,  which  can  be  much  larger  than  that  due  to  the 
boundary  layer  itself. 

Very  roughly,  drag  due  to  flow  separation  may  be  repre- 
sented in  a  two-dimensional  velocity  profile  with  the  reference 
frame  fixed  on  the  ship,  as  shown  in  Figure  6. 


+00  +6 

pudJ^  -  u)dy  - 


pu(U   -  u)dy  (10a) 


From  Figure  3,  it  can  be  seen  that  for  a  separated  wake,  the 
contribution  to  this  integral  from  separation  along  is  approx- 
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This  quantity  may  be  quite  large.   If  flow  separation 
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occurs,  then  the  form  drag  will  increase  by  this  large  amount 
and  may  increase  to  a  value  many  times  that  of  the  form  drag 
due  only  to  an  attached  boundary  layer. 

An  additional  component  of  drag  which  is  related  to  the 
form  of  the  hull  is  the  difference  in  frictional  resistance 
between  the  wetted  surface  of  the  ship  in  a  plane  uniform  flow 
and  that  in  the  curved  form  of  the  hull.   Because  flow  is 
forced  around  the  ship  in  a  curved  path,  the  pressure  distri- 
bution causes  an  increase  in  the  velocity  near  midships  and  a 
decrease  near  the  ends.   The  increase  in  velocity  around  the 
midship  region  has  a  greater  effect  on  the  frictional  drag 
than  the  reduction  in  velocity  at  the  ends.   The  result  is  an 
overall  increase  in  frictional  resistance  over  that  of  plane 
uniform  flow.   In  addition  to  the  velocity  field  changes,  the 
path  is  longer  along  a  curved  hull  than  along  a  flat  plate  of 
the  same  area.   These  effects  are  both  more  pronounced  for 
shapes  with  high  volumetric  coefficients  than  for  very  slender 
forms  with  low  volumetric  coefficients.   Since  this  friction 
drag  is  related  to  the  form  of  the  hull,  it  is  included  in 
form  drag. 

If  form  drag  were  plotted  as  a  function  of  volumetric  co- 
efficient and  then  extrapolated  to  zero  volumetric  coefficient, 
the  value  of  the  drag  at  the  resulting  intercept  point,  al- 
though included  in  form  drag  by  conventional  test  procedures, 
cannot  be  due  to  the  amount  of  displacement  contained  in  the 
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hull  form  because  the  hull  can  have  no  thickness  at  a  zero 
value  of  volumetric  coefficient.   Since  frictional  resistance 
is  not  measured  directly  in  the  towing  tank,  but  is  calculated 
by  whatever  frictional  drag  formula  is  adopted  by  the  test 
procedure,  then  the  value  of  the  form  drag  at  zero  volumetric 
coefficient  constitutes  primarily  the  difference  between  the 
value  predicted  by  the  adopted  friction  formula  and  that 
actually  contributing  to  total  drag.   Where  they  are  signifi- 
cant, any  other  errors  in  the  measurement,  calculation,  and 
scaling  of  frictional  drag  would  also  end  up  in  this  category. 


-30- 


IV.   THE  INTERACTION  BETWEEN 
WAVEMAKING  AND  FORM  DRAG 


The  thickness  of  the  boundary  layer  is  dependent  on  the 
pressure  gradient  along  the  path  of  the  flow.   The  pressure  in 
general  increases  up  to  the  stem,  decreases  toward  midships, 
increases  up  to  the  stern  and  decreases  back  to  ambient  in  the 
wake.   Flow  is  thereby  decelerated  ahead  of  the  ship,  acceler- 
ated from  the  stem  to  midships,  decelerated  from  midships  aft 
and  accelerated  slightly  into  the  wake.   The  fluid  entrained 
in  the  boundary  layer  is  then  flowing  against  the  pressure 
gradient  from  the  midship  section  aft.   The  boundary  layer 
will  therefore  thicken  at  a  greater  rate  aft  than  it  will  in 
the  accelerating  flow  forward.   The  pressure  decrease  from  the 
stem  to  midships  is  greater  than  the  pressure  recovery  from 
midships  to  the  stern.   This  causes  the  viscous  pressure  drag 
mentioned  in  the  previous  section.   It  has  the  somewhat  re- 
deeming virtue,  however,  that  the  adverse  pressure  gradient 
against  which  the  fluid  in  the  boundary  layer  must  flow  is 
decreased  in  a  gentle  and  continuous  manner  from  that  which 
would  exist  with  a  complete  recovery  of  the  pressure  at  the 
stern. 

All  of  the  foregoing  refers  only  to  the  pressure  distri- 
bution which  would  exist  without  waves.   When  the  pressure 
distribution  of  a  wave  system  is  superimposed  upon  the  pressure 
distribution  discussed  in  the  previous  paragraph,  skepticism 
arises  concerning  the  predictability  of  the  interaction  between 
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two  completely  different  phenomena  which  are  each  in  their  own 
right  relatively  unpredictable. 

One  of  the  most  insidious  interactions  lurking  in  the 
towing  tank  might  be  the  initiation  of  flow  separation  by  the 
pressure  field  of  the  wave  system.   If  the  flow  were  on  the 
verge  of  separation  near  the  stern  of  a  model  and  if  the  model 
were  then  run  at  a  speed  corresponding  to  a  hollow  in  the  re- 
sistance curve,  the  crest  of  a  transverse  wave  would  be  posi- 
tioned at  the  stern,  adding  a  pressure  peak  near  the  end  of  an 
already  adverse  pressure  gradient,  whereupon  flow  separation 
might  occur.   This  sudden  increase  in  drag  might  then  fill  in 
all  or  part  of  an  anticipated  hollow  in  the  resistance  curve. 
If  the  speed  were  then  changed  to  that  corresponding  to  a  hump 
in  the  anticipated  resistance  curve,  a  transverse  wave  trough 
would  be  positioned  near  the  stern,  reducing  the  pressure  and 
accelerating  the  flow  toward  the  stern.   Any  region  of  previ- 
ously separated  flow  might  be  partially  or  completely  re- 
attached, suddenly  reducing  the  drag  and  removing  all  or  part 
of  the  expected  hump  in  the  drag  curve. 

In  addition  to  the  effect  of  the  wave  system  on  flow 
separation,  flow  separation  may  in  turn  affect  the  wave  system. 
Mr.  P.N.  Joubert,  in  his  comments  on  reference  9,  cites  an 
example  of  what  appeared  to  be  a  wave  induced  by  flow  separa- 
tion on  the  12  meter  yacht  "Valiant". 

In  this  way,  the  resistance  curve  of  a  model  test  might 
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be  drastically  altered  by  an  interaction  effect.   If  the  model 
were  then  scaled  up  to  a  full  size  ship,  the  very  large  Rey- 
nold's number  might  prevent  flow  separation  even  in  the  pres- 
ence of  the  same  wave  profile  as  was  modelled.   The  resistance 
of  the  full  size  vessel  might  then  be  entirely  different  from 
its  model  test. 

Even  without  flow  separation,  similar  effects  might  be 
produced  to  a  lesser  degree  simply  by  the  influence  of  the 
wave  system  pressure  field  on  the  thickening  and  thinning  of 
the  boundary  layer  over  the  afterbody.   R.L.  Townsin  [25]  said 
that  the  reduction  in  the  non-dimensional  boundary  layer  thick- 
ness, 6/L,  may  be  as  great  between  an  8-foot  model  and  a  28 
foot  model  as  it  is  between  the  28-foot  model  and  the  full 
size  ship.   It  is  also  agreed  in  the  introduction  to  reference 
25  that  separation  scale  effects,  if  they  exist,  may  appear 
over  the  range  of  model  size. 
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V.   THE  DATA  BASE 

The  selection  of  a  methodical  series  of  ship  model  tests 
for  use  as  a  data  base  for  establishing  trends  and  functional 
relationships  between  wave  drag  and  other  parameters  is  limited 
by  several  considerations.   The  first  of  these  considerations 
is  the  absolute  size  of  the  model. 

Regardless  of  the  causes,  the  average  magnitude  of  the 
errors  encountered  in  tank  test  results  increases  strongly 
with  decreasing  model  size  and  Reynold's  number.   Figure  7 
shows  this  relationship.   It  is  suspected  by  the  author  that 
irregularities  in  tank  test  data  may  result  from  the  interac- 
tions between  the  wave  system  pressure  field  and  the  boundary 
layer  along  the  afterbody,  as  discussed  in  the  previous  sec- 
tion.  Uniformity  in  the  behavior  of  the  boundary  layer  is 
improved  by  increasing  the  Reynold's  number  and  since  the 
speed  and  kinematic  viscosity  are  not  available  as  control 
variables  for  this  purpose,  then  increasing  the  size  of  the 
model  will  produce  the  same  result.   To  insure  that  wave 
effects  are  measured  in  a  flow  regime  similar  to  that  of  a 
full  size  ship,  it  is  extremely  important  that  any  irregular- 
ities due  to  small  model  size  be  avoided.   For  this  reason, 
all  ship  model  test  data  for  models  less  than  10  feet  long 
were  rejected  as  unreliable. 

When  one  requires  the  size  of  the  model  to  be  larger  than 
10  feet,  testing  becomes  extremely  expensive  and  methodical 
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series  for  large  models  in  recent  years  tend  to  include  only 
small  numbers  of  models.   This  brings  us  to  the  next  consider- 
ation, that  of  statistical  significance.   Data  populations 
with  only  a  few  variations  in  some  parameters  of  interest  may 
lead  to  results  which,  regardless  of  their  correctness,  will 
not  stand  up  under  any  statistical  measure  of  credibility.   An 
example  of  this  problem  is  provided  by  the  BSRA  Trawler  Series 
(Part  I)  [17],  wherein  a  total  of  only  nine  models  were  built. 
Four  models  with  varying  slenderness  at  a  constant  beam  to 
draft  ratio  were  tested  along  with  five  models  of  varying  beam 
to  draft  ratio  at  a  constant  slenderness.   With  no  variation 
in  prismatic  coefficient  and  so  few  in  other  parameters,  Mr. 
G.S.  Selman  commented,  "It  surely  needs  a  lot  of  imagination 
to  produce  from  four  models  such  wonderful  curves,  as  shown.." 
[18].   Unless  criticism  of  this  nature  is  to  be  invited,  a 
larger  data  base  is  required.   This  methodical  series  could 
not  be  used,  therefore,  even  with  a  model  size  of  twelve  and 
a  half  feet  and  a  scale  ratio  of  only  one  1/12. 

When  other  parameters  are  desired,  the  BSRA  Trawler  Series 
becomes  more  attractive.   Although  the  number  of  data  points 
in  any  one  parameter  is  limited,  changes  in  LCB,  C_.  and  load 

a 

condition  are  available  if  Part  I  and  Part  II  [17  and  18]  are 
both  used.   Valuable  insights  have  been  gained  from  multiple 
regression  analysis  of  the  entire  series  [17]. 

Series  60  [24]  was  investigated  and  even  with  large  numbers 
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of  large  models,  could  not  be  used  because  the  range  of  pris- 
matic coefficient  did  not  include  values  sufficiently  low  to 
be  applied  to  many  vessels  other  than  merchant  ships.   In 
addition,  the  parent  lines  of  Series  60  possess  shoulders 
which,  although  appropriate  for  merchant  ships,  might  intro- 
duce the  additional  complexities  of  the  fore  and  aft  shoulder 
wave  systems  similar  to  those  of  the  double-wedge-shaped  body 
investigated  by  Wigley  [27] . 

After  consideration  of  several  other  methodical  series, 
the  Taylor  Standard  Series  [7  and  22]  was  selected  for  a  number 
of  reasons.   The  models  are  all  about  20  feet  long  and  average 
about  a  ton  in  displacement.   The  number  of  data  points  avail- 
able in  the  variation  of  volumetric  coefficient,  prismatic 
coefficient  and  speed  length  ratio  are  sufficiently  large  to 
clearly  represent  functional  relationships.   The  lines  of  the 
parent  model  for  the  Taylor  Standard  Series  are  presented  in 
Figure  8.   Both  LCB  location  and  midship  section  coefficient 
are  held  constant  in  the  Taylor  Standard  Series. 

An  additional  attraction  to  the  Taylor  Series  as  a  data 
base  is  the  manner  in  which  the  sectional  area  curve  and  the 
lines  of  the  parent  are  changed  to  produce  the  desired  varia- 
tions in  form  coefficients.   The  sectional  area  curve  is 
defined  by  a  fifth  order  polynomial  which  is  related  directly 
to  the  prismatic  coefficient  and  which  is  proportioned  for  the 
variations  in  the  series  in  a  smooth  reproducible  manner.   The 
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lines  are  changed  from  the  parent  graphically  in  a  process 
consistent  with  the  changes  in  the  sectional  area  curve  and 
produce  very  fair  hull  forms  without  shoulders  or  other  areas 
of  high  local  curvature. 

Almost  200  models  were  tested  from  1902  to  1914,  of  which 
120  were  published  in  The  Speed  and  Power  of  Ships  [22] .   Al- 
most 160  of  these  tests  were  later  faired  into  the  data  pre- 
sented in  the  Gertler  re-analysis  of  the  original  test  data 
[7].   Although  the  data  is  old,  it  is  not  out  of  date  for  the 
purposes  of  this  thesis.   In  fact,  a  single  methodical  series 
of  the  magnitude  of  the  Taylor  Standard  Series  will  probably 
never  again  be  attempted. 

The  data  published  in  the  three  editions  of  Taylor's  The 
Speed  and  Power  of  Ships  [22]  included  only  beam  to  draft 
ratios  of  2.25  (Series  21)  and  3.75  (Series  22)  between  which 
interpolation  was  recommended.   Series  20  contained  the  data 
for  the  intermediate  value  of  beam  to  draft  ratio  of  2.92. 
Although  Series  20  was  completed  prior  to  Series  21  and  22,  it 
was,  for  some  reason  unknown  to  this  author,  not  published 
along  with  Series  21  and  22  as  an  intermediate  value  from  which 
to  interpolate  beam  to  draft  ratio  effects.   In  addition,  Series 
20  was  not  extended  in  1913  and  1914  as  were  Series  21  and  22. 
Taylor  did  not  explain  why  this  data  from  some  38  models  was 
left  out,  but  only  mentions  in  passing  that  it  confirmed  to 
some  extent  the  validity  of  linear  interpolation  between  Series 
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21  and  22. 

"A  Reanalysis  of  the  Original  Test  Data  for  the  Taylor 
Standard  Series"  [7]  was  published  in  1954  by  Morton  Gertler. 
Series  20  was  included  in  these  data  and  was  used  to  fair  the 
data  and  provide  an  intermediate  value  of  beam  to  draft  ratio 
of  3.0.   The  original  Taylor  data  was  corrected  for  many  ef- 
fects and  updated  for  modern  conventions  in  this  reanalysis. 
For  this  reason,  the  data  presented  in  Gertler* s  reanalysis 
were  used  in  this  thesis  with  only  one  reservation.   The 
behavior  of  the  curves  of  residuary  drag  coefficient  with  that 
value  at  V//L  =  .5  subtracted  out  and  then  multiplied  by  the 
wetted  surface  plotted  against  volumetric  coefficient  for  each 
beam  to  draft  ratio,  speed  length  ratio  and  prismatic  coeffi- 
cient reveal  a  very  unusual  behavior  for  the  curve  representing 
the  beam  to  draft  ratio  of  3.0.   While  the  curves  for  the  B/T 
of  2.25  and  3.75  behave  similarly,  the  curve  for  the  B/T  of 
3.0  usually  starts  out  shooting  above  the  other  two  and  then 
crossing  over  the  other  two  at  a  higher  value  of  volumetric 
coefficient.   This  strange  behavior  of  the  data  for  B/T  of  3.0 
is  evident  at  many  values  of  prismatic  coefficient  and  speed 
length  ratio  when  plotted  carefully  on  a  large  scale.  Examples 
are  plotted  in  Figures  9,  10  and  11.   Because  of  this  unusual 
and  unexplainable  behavior,  the  data  for  B/T  =  3.0  is  not  used 
in  this  thesis.   It  is  possible  that  a  difference  in  the  be- 
havior between  the  data  for  Series  20  and  that  of  Series  21 
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and  22  may  reflect  some  important  problem  with  the  models  or 
the  tests  that  was  suspected  or  recognized  by  Admiral  Taylor 
himself  and  was  for  that  reason  omitted  from  his  publications 


-44- 
VI.   WAVE  DRAG  IN  THE  DATA  BASE 

Since  the  Gertler  reanalysis  of  the  Taylor  Standard 
Series  was  chosen  as  a  data  base,  the  range  of  length  to  draft 
ratio  for  which  data  was  obtained  will  be  examined  to  estimate 
the  magnitude  of  the  effect  of  the  distribution  of  displace- 
ment in  depth  for  reasonable  wave  numbers  and  speed  length 
ratios. 
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BT    V    1      Cv 


L2    L3  C  C     C  C 
p  m    p  m 


(12) 


L*  =  (SESE)BT  (13) 

Lv 


<?>  ■  4^]  $  (i5) 


'I'max  ^--rmr1  (3-75)  "  49-2  (16) 


(£)  .   =a/(-4)8>'-9,23)(2.25)  =  12.9  (17) 

T  mm   V   (.006) 


The  variation  in  the  velocity  potential  decay  rate  seen  by 
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this  data  base  is  then 


*<<!>max>=V1/(F)<L/T)=    V1/49-2F2 


♦  ((&.)-♦    e-l/(P2)(VT)    .    #    e-l/12.9F2 


Since  speed  length  ratios  from  0.5  to  1.2  were  included 
in  the  data  base,  the  variation  of  the  effects  of  the  length 
to  beam  ratios  included  in  the  data  base  on  the  wave  drag  will 
be  estimated.   Froude  number  is  proportional  to  the  speed 
length  ratio  as  follows: 


_     V     1.68  9VK    ,1-689, VK     9Q7fiVK  ,9m 


At  a  VK//L  =  1.2,  Froude  number  is  .357.   Then, 


tuh         ,    -5—  1.2)  =  0   e-1/(-357)2(49-2»  =  .853*   (21) 
T   T  max'  /L  yo  Yo 


a//L»       V    .  0,    .   -1/  .357  2  12.9      ....    ,00« 
4>  (  (™)  •  /  ~^=  =  1.2)  =  *  e  =  .544*   (22) 

f   T  mm  fL  Yo  *o 


It  may  be  seen  that  at  this  speed  length  ratio,  there  is 
very  little  variation  over  depth  of  the  velocity  potential  for 
the  larger  length  to  draft  ratio.   For  the  smaller  length  to 
draft  ratio,  however,  the  velocity  potential  is  almost  reduced 
in  half  at  the  draft  of  the  ship.   A  given  redistribution  in 
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depth  of  sectional  area  or  displacement  would  then  presumably 
cause  a  larger  change  in  wave  drag  for  the  smaller  length  to 
beam  ratio  than  for  the  larger. 

At  a  speed  length  ratio  of  1.0,  the  Froude  number  is: 


F  =  (.2976)-£  =  .2976  (23) 


Then, 


w/L^  v  i    n\  a      -1/    .2976    2    49.2  _QC.       ,„,, 

<f>  (  (m)  /    -/¥=  =    1.0)    =    <b   e  =    .795<t      (24) 

yuT  max      /L  o  Yo 


"fWr1'«»  =  V1/u"6,!,12'"-"'V25' 


Again  it  may  be  seen  that  the  variation  of  4>(z)  is  quite 
large  for  the  smaller  length  to  draft  ratio.   The  variation  of 
4>(z)  for  the  larger  length  to  draft  ratio,  although  greater  at 
this  speed,  is  still  not  a  large  fraction  of  <j)  . 

At  a  speed  length  ratio  of  . 8 , 


VK 
F  =  (.2976)-t£  =  .2381  (26) 


and 


f((L)    ,  V  e-l/(.2381)2(49.2)  _ 

y   T  max'  /L  vo  Yo 
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A//L,      V     Qv    ,   -l/(. 2381) 2 (12.9)     occ,    ,00* 
^min'  7r=  '8)  =  V  =  '255*o   (28) 

At  this  speed  length  ratio,  the  variation  of  <\>(z)    is 
large  even  for  the  large  length  to  draft  ratio.   At  the  same 
time,  however,  wave  drag  has  decreased  to  a  small  value  such 
that  the  absolute  change  in  wave  drag  due  to  a  given  redistri- 
bution in  displacement  in  depth  would  probably  not  be  very 
large. 
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VII.   HYPOTHESIS 

It  is  hypothesized  that  the  residuary  drag  may  be  decom- 
posed into  three  components:   wave  drag,  form  drag  and  skin 
friction  error. 


RR  "  R0  +  RI  +  Rw  (29» 


where    R  =  wave  drag  =  -^-pV2L2C 
w  3    2      w 

C   =  wave  drag  coefficient  based  on  waterline  length 

R,  =  form  drag  =  ypV2L2C, 

C,  =  form  drag  coefficient  based  on  length 

R  =  skin  friction  error  =  -=-pV2L2C 
o  2K     o 

C  =  skin  friction  error  coefficient  based  on  length 
o 


It  is  further  hypothesized  that  the  wave  drag  will  be  a 
power  law  function  of  the  displacement,  V*,  of  the  ship  weighted 
by  its  distribution  in  depth  and  an  unknown  function,  F_,  of 
the  longitudinal  prismatic  coefficient  of  the  weighted  dis- 
placement and  the  speed  length  ratio. 


\  =  !?V2I-2CW  (30) 


cw=  K^'VV'  7E'1  (31) 


where    n  =  exponent  to  be  determined 
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F..  =  unknown  function  to  be  determined 


Because  it  was  expected  that  n  would  be  approximately  2.0, 
a  preliminary  function  was  defined  such  that: 


c„=  "P'2"f2(<.  >  (32) 


F  could  then  be  more  easily  determined  and  F3  could  be  related 
back  to  F?  when  n  became  known  more  precisely. 


c  =  F  (V\n  =  F  (V*  2  (33) 

w     3"'      2VL3' 


So, 


F   =  P  (V\(2-n)  (34) 

*3     2VL3' 


In  addition,  it  is  expected  that  the  form  drag  will  be 
proportional  to  the  real  displacement,  V,  and  also  proportional 
to  some  unknown  function,  F^of  the  beam  to  draft  ratio  and  the 
real  (unweighted)  longitudinal  prismatic  coefficient. 


Rx  =  Ipv'l2^  (35> 


C      =    rXl  TF    (-     CI  <36> 

Cl         1LtJ  ltl  V      pJ 
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The  skin  friction  error  is  defined  as  the  value  of  the 
form  drag  when  it  is  extrapolated  to  zero  volumetric  coeffi- 
cient.  It  is  defined  in  such  a  way  that  it  represents  that 
part  of  the  form  drag  which  cannot  depend  on  the  volumetric 
coefficient  of  the  hull  form.   It  would  seem  at  first  glance 
that  since  the  hull  form  can  have  no  volume  to  shape,  the  skin 
friction  error  should  not  depend  on  any  form  coefficient  what- 
soever.  The  other  form  coefficients  are  not,  however,  forced 
by  the  definition  of  skin  friction  error  to  be  zero  and  it  is 
expected  that  some  of  the  errors  included  within  the  definition 
may  vary  with  beam  to  draft  ratio  and  prismatic  coefficient. 
In  addition,  since  the  real  skin  friction  is  a  function  of  the 
hull  form,  the  skin  friction  error  will  in  general  also  be  a 
function  of  hull  form  as  follows: 


R  =  ipV2L2C  (37) 


o    2 K     o 


co  =  Vf<  V  (33> 
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VIII.   DEFINITION  OF  NEW  VARIABLES 

Because  the  variation  of  the  velocity  potential  of  the 
wave  system  over  depth  depends  upon  both  the  Froude  number  and 
the  depth  itself,  new  variables  must  be  defined  which  relate 
this  dependency  to  the  speed  and  form  of  the  ship.   In  order 
to  include  the  effect  of  the  vertical  distribution  of  dis- 
placement on  the  velocity  potential  with  ship  speed,  a  new 
displacement  parameter  is  defined  which  is  an  integral  over 
draft  with  the  integrand  weighted  by  depth  and  Froude  number. 

The  area  of  each  section  may  be  weighted  as  follows: 

z=T 
A*x)  =    y(z,x)e~gz/V  dz  (39) 

z=0 

The  weighted  displacement  is  then: 


V*  = 


x=L 

t 


x= 

=L    z= 

=  T 

* 

' 

— 

[ 

J 

x= 

=0    z= 

=  0 

A(x)dx  (40) 

=  0 


(z,x)e"g/V  zdz]dx  (41) 


where  y(z,x)  are  the  offsets  determined  graphically  from  ref- 
erence 7 ,  and  e  g      represents  the  weighting  factor  for  the 
contribution  of  the  local  displaced  volume  to  the  disturbance 
of  the  velocity  potential. 

It  may  be  seen  that,  in  general, 
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V   <  V  (42) 


and  that  for  high  speeds,  V*  approaches  the  real  displacement. 
As  discussed  in  the  section  on  wave  drag  and  the  data  base, 
all  of  the  displaced  volume  contributes  approximately  the  same 
toward  the  disturbance  of  the  free  surface  at  very  high  speeds 

In  a  manner  similar  to  the  derivation  of  standard  coeffi- 
cients of  form,  a  new  set  of  weighted  coefficients  may  be 
derived  which  include  in  them  their  relative  contributions  to 
the  disturbance  of  the  velocity  potential. 

Define  a  weighted  midship  section  area, 


z=T 

* 

A  = 
m 

z=0 


y(z)e'gz/V  dz  =  a|x)  at  x  =  51 


(43) 


Then  a  weighted  longitudinal  prismatic  coefficient  may  be 
defined: 


*    y* 
C   =  — —  (44) 

P    A™L 

m 


The  weighted  volumetric  coefficient  would  be  written, 


*    v* 

cv  "  b  (45) 


The  longitudinal  center  of  weighted  buoyancy  is: 
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x=L 


x=L 


x=L 


LCB 


x  A(x)dx 


A(x)dx  =  775- 


1^ 

V* 


x  A(x) dx 


(46) 


x=0 


x=0 


x=0 


Although  a  variation  in  LCB  was  not  desired  and  LCB  is 
fixed  amidships  in  the  Taylor  Standard  Series,  LCB*  will,  in 
general,  change  with  speed.   The  Taylor  Series  ships  tend  to 
have  U-shaped  sections  forward  and  V-shaped  sections  aft.   At 
low  speeds,  more  V*  will  be  lost  forward  than  aft  and  will 
result  in  a  migration  of  LCB*  aft  at  low  speeds.   To  insure 
that  this  variation  remained  small  enough  to  ignore,  LCB*  was 
defined  and  calculated. 
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IX.   DATA  REDUCTION 

Within  the  data  base  of  the  Taylor  Standard  Series,  a 
representative  subset  was  chosen  with  enough  data  points  to 
clearly  establish  trends  between  parameters.   The  subset  of 
ships  contains  70  ships  with  the  following  range  of  variation 
of  coefficients  of  form: 


|  :  2.25,  3.75 


CV:  .001,  .002,  .003,  .004,  .005,  .006,  .007 


C  :  .48,  .54,  .60,  .65,  .70 
p 


Fractional  offsets  were  determined  graphically  for  each  vari- 
ation of  C   for  10  stations  and  the  0,  .1,  .2,  .3,  .4,  .5,  .6, 
.8  and  1.0  waterlines. 

For  those  residuary  drag  curves  in  reference  7  which  are 
flat  at  the  speed  length  ratio  of  .5,  the  residuary  drag  co- 
efficient at  this  speed,  CD(.5),  is  assumed  to  be  essentially 
form  drag.   The  form  drag  coefficient,  C, ,  based  on  length 
squared,  may  then  be  written  as: 


Cl  E  1/2RPV^  (47) 


In  reference  7,  the  wetted  surface  coefficient  is  defined  as: 
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C0  =  "7§r  (48) 


's  "  7vn 


and 


cv E  h  (49) 


so  that  the  combination  of  coefficients 


. Rr(.5)     q     y      RR(.5) 

cR(.5)cs/c;  =  (1/2  pVis)  (^  (  £r>  =  iTf^v^  ■  Cl  (50) 


is  equal  to  the  form  drag  coefficient,  C, . 

The  form  drag  coefficient,  C. ,  was  then  plotted  against 

C   for  various  values  of  C   and  B/T.   The  form  drag  coefficient 
v  P 

was  found  to  be  essentially  linear  in  C   and  was  rewritten  in 

J  v 

the  form 


Ci  =  FlCv  (51) 


where  F,  is  simply  the  slope  of  the  curve  of  form  drag  coeffi- 
cient as  a  function  of  C  .   F,  is,  however,  an  unknown  function 

v     1 

of  C   and  B/T.   The  form  drag  is  then 


Rx  =  |pV2L2CvF1(Cp,  |)  (52) 


The  intercept  of  the  curve  of  form  drag  coefficient  versus 
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C  extrapolated  to  C  =0  is  taken  as  the  skin  friction  error 
v  v 

coefficient  based  on  length.   So, 


F   =  CD  Ce/C~  =  =  yo  „2t  2    =   C  (53) 

o     RQ  S   v    1/2  pV  L      o 


where  F   is  not  a  function  of  C  ,  but  is  an  unknown  function 

of  C  and  B/T. 
P 

The  skin  friction  error  is  then, 


R0  =  |pV2L2F0(Cp,  |)  (54) 


Examples  of  the  plots  of  form  drag  coefficient  versus  C   are 
provided  in  Figures  12,  13  and  14.   Examples  of  plots  of  the 
form  drag  coefficient  versus  C   are  provided  in  Figures  15, 
16  and  17.   Although  these  plots  include  data  for  the  beam  to 
draft  ratio  of  3.0,  these  data  were  not  used  in  the  derivation 
of  form  drag  functions. 

Linear  regression  analyses  were  done  graphically  to  de- 
termine reasonable  fits  for  the  functions  F, (C  ,B/T)  and 

1      p 

F    (C    ,    B/T).      The    following    relationships   were   chosen: 


F      =    .0027    +    (C      -    .5)  [.021+   (|  -    2.25)  (.02033)]  (55) 

op  T 


F,    =    14    +    (|  -    2.25)  [4    +    20.7(C      -    .5)]  (56) 

IT  p 


-57- 


Plots  of  F   versus  C   are  provided  for  three  values  of 

0  P 

B/T  in  Figure  18. 

Plots  of  F,  versus  C   are  also  provided  for  three  values 

1  p 

of  B/T  in  Figure  19  with  some  data  points  included. 
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Plot   of    Form     drag  as    rt     function    of 
Volumetric     coefficient    for.     various 
Values     of     Prismatic     Coefficient 
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Plot  of    Form  drag    coefficient 
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A  computer  program  was  written  to  do  the  necessary  calcu- 
lations and  is  presented  in  one  version  in  Appendix  A.  Offsets 
were  read  and  interpolated  to  provide  offsets  at  .7  and  .9 
waterlines  which  were  not  available  in  reference  7.   Complete 
full  size  offsets  were  generated  for  ten  stations  and  ten 
waterlines  for  each  of  the  70  ships.   From  these  offsets,  V, 

S,  C  ,  C  ,  C   and  C„  were  calculated.   All  integrals  were  done 
v   p   m      S  3 

using  Simpson's  rule.   Calculations  were  done  for  vessels  of 
190  feet  in  length. 

For  each  of  the  70  ships  the  weighted  form  coefficients 
were  calculated  for  15  speeds,  providing  1050  weighted  forms 
from  the  1050  ship-speed  combinations.   Although  this  thesis 
contains  formulae  for  predicting  form  drag,  they  were  pro- 
grammed only  as  a  means  of  comparing  predicted  form  drag  to 
CD(V//L  =  .5).   To  prevent  the  validity  of  one  hypothesis 
from  depending  on  the  validity  of  another  hypothesis,  the 
values  of  CD(V//L  =:  .5)  punched  from  the  Taylor  Series  were 
used  as  form  drag  in  the  calculation  of  wave  drag.   F~  was 
then  calculated  as 


F2  "  tCR"  CR(-5)1(£)(C02  (57) 


The  product  of  the  first  computer  program  was  both  printed 
and  punched  for  further  use.   Printed  outputs  for  a  number  of 
example  ships  are  included  in  Appendix  C.   For  each  of  the 
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1050  ship-speed  combinations,  the  following  vital  parameters 
were  punched  on  cards  for  later  analysis: 


NAME 
SLR 
CCR 
CR5 

CVST 
CVT 

CPS1 

CPN 

SL2 

BDR 

F2 


QUANTITY 
V//L 
CR(V//L) 
CR(.5) 

C*  x  10 3 

v 

C   x  103 
v 

CP* 

CP 

S/L2 

B/T 

F„ 


The  printed  output  contains  additional  information,  including 
Fo'  Fl'  S/  V*'  Froude  number,  LCB*,  C*  and  C*. 
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X.   THE  REGRESSION  ANALYSIS 

The  Statistical  Package  for  the  Social  Sciences  [2]  was 
used  to  conduct  all  regression  analyses  and  correlation  studies. 
This  very  powerful  analytical  tool  provides  extremely  easy 
interfacing  and  an  output  containing  every  imaginable  measure 
of  statistical  significance.   Simple  regression,  multiple  re- 
gression and  stepwise  inclusion  multiple  regression  processes 

are  available. 

■ 

The  1050  ship-speed  combinations  were  filtered  to  elimi- 
nate misleading  data  and  shakey  extrapolations.   A  ship  with 
Cv  =  .007  and  CD  =  .7 ,  for  example,  at  a  V//L  =  1.2  will  have 
a  residuary  drag  coefficient  so  large  that  it  would  plot  off 
the  page  in  reference  7  if  the  extrapolation  were  attempted. 
Over  100  ship-speed  combinations,  although  extrapolated  as 
data  in  reference  6,  were  eliminated  for  these  reasons. 

Correlation  studies  were  conducted  in  an  attempt  to  deter- 
mine n  in  the  equation 


Cw  ■  F3(Cv»n  "  F2(Cv>2  (58) 


for  the  entire  1050  ship-speed  combinations.   After  investiga- 
ting various  values  of  n  by  computer  and  calculating  best  values 
of  n  by  hand,  n  was  fixed  at  a  value  of  1.8  for  all  speeds,  so 
that: 
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C*"*3lCVU*  <59) 


It  should  be  mentioned  at  this  point  that  a  value  of  n  =  1.8 
very  slightly  favors  high  speeds  from  V//L  =  1.0  to  1.2.   For 
lower  speeds,  a  somewhat  lower  value  of  n  would  have  provided 
a  better  fit.   For  the  purposes  of  this  thesis,  however,  a 
constant  value  of  n  =  1.8  was  used  for  all  speeds. 

Due  to  the  irregular  manner  with  which  the  relationships 
between  the  form  coefficients  and  residuary  drag  change  with 
speed  length  ratio,  a  separate  multiple  regression  analysis 
was  run  for  each  of  the  14  speed  length  ratios  from  .55  to  1.2. 
In  addition  to  the  variable,  CP*,  a  series  of  five  sine  func- 
tions of  CP*  was  defined  and  made  available  for  stepwise 
inclusion  into  each  multiple  regression  analysis  as  necessary 
to  obtain  the  desired  fit.   Since  the  values  of  F~  were  punched 
in  the  data  rather  than  F-^,  it  was  necessary  to  compute  F~  as 
a  variable  within  the  regression  analysis  program  for  every 
sample  along  with  the  terms  in  the  Fourier  sine  series.   Ap- 
pendix B.l  contains  selected  listings  and  output  printouts  for 
one  of  the  14  regression  analyses  runs.   Variables  were  defined 
as  follows: 

CPS1  =  CP* 

CVST  =  CV*  x  10 3 

CPSV  E  (CPS1  -  .48)/. 27 


-70- 


CPSA  =  sin(TTCPSV) 

CPSB  E  sin(27TCPSV) 

CPSC  E  sin(37rCPSV) 

CPSD  =  sin(47TCPSV) 

CPSE  =  sin(5TrCPSV) 

F3  E  F2(CVST) *2 

The  Fourier  sine  series  was  limited  to  five  terms  to 
avoid  the  representation  of  data  scatter  in  the  final  function 
Each  stepwise  inclusion  multiple  regression  analysis  then  re- 
gressed F^  with  CPS1  and  the  five  sine  terms.   A  constant  is 
automatically  added  into  the  regression  as  the  zeroth  order 
term.   The  final  form  of  the  universal  function,  F-.,  is  then: 


F-  =  CK  +  C,CPS1  +  C.CPSA  +  C0CPSB  +  C^CPSC 
3  1         A         B         C 


+  C^CPSD  +  C^CPSE  (60 

D         hi 


The  resulting  values  of  the  constants  are  listed  in  Ap- 
pendix B.2  along  with  various  measures  of  confidence  and 
statistical  fit.   An  example  of  a  typical  listing  of  SPSS  con- 
trol cards  used  is  also  provided  in  Appendix  B.l  along  with 
parts  of  the  SPSS  output  for  that  particular  run.   It  should 
be  noted  that  the  F-.  in  the  regression  analyses  computer  pro- 
gram is  larger  than  the  F^  defined  in  previous  sections  by  a 
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2 

factor  of  (1000)*   =  3.98  due  to  the  definition  of  CVST, 


CVST  =  CV*  x  10 3  (61) 


To  use  this  universal  function  for  the  prediction  of  wave  drag, 
F-.  must  be  divided  by  (1000)  *2  or: 


F3    =    F2(CV*)*2    =    F2(CVST)*2/(1000)  '2  (62 


2 

[Tables   of   F^   =   F~ (CV*) *      are   provided    in   Appendix   B.3. 
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XI.   WAVE  DRAG  CONCLUSIONS 

For  each  value  of  V//L,  F.,  constitutes  an  estimate  of  a 
universal  wave  drag  function  of  C*  alone  and  represents  wave 
drag  coefficient  per  unit  weighted  volumetric  coefficient 
raised  to  the  1.8  power.   With  the  use  of  the  weighted  coef- 
ficients of  form,  C*  and  C*,  and  with  the  implementation  of 
the  universal  wave  function,  F-.,  the  representation  of  wave 
drag  is  simplified  considerably.   The  entire  Taylor  Standard 
Series,  for  example,  is  collapsed  onto  one  sheet  of  paper  as 
represented  in  Figure  20.   A  family  of  curves  of  F_  versus  C* 
is  plotted  in  this  illustration  with  one  curve  for  each  value 

of  V//L.   It  is  reassuring  that  the  curves  of  F_,  versus  C*  for 

3         p 

various  values  of  V//L  show  minima  similar  to  those  found 
using  the  standard  coefficients  of  form.   Figure  21  illustrates 
these  minima  for  the  Taylor  Standard  Series  data  in  reference 
7. 

The  value  of  the  exponent  of  displacement  (or  C  )  would 
be  expected  to  be  2.0  from  Mitchell's  [14]  linear  ship  theory. 
Weinblum  [26]  states  that  wave  drag  was  approximately  propor- 
tional to  displacement  to  the  1.6  power.   With  the  form  drag 
subtracted  out  of  residuary  drag,  it  is  concluded  from  the 
results  of  this  thesis  that  wave  drag  is  approximately  propor- 
tional to  the  displacement  raised  to  the  1.8  power.   Thus  half 
of  the  discrepancy  between  the  theoretical  value  of  2.0  and 
the  previously  observed  value  of  1.6  has  been  explained  and 
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removed . 

Figure    22    illustrates    the   effect   of   both   C      and   of   C      on 

p         v 

the  residuary  resistance.   Some  sample  curves  of  wave  drag 
coefficient  predicted  by  the  universal  wave  drag  function  and 
by  the  Taylor  Standard  Series  are  provided  for  comparison  in 
Appendix  C.   From  these  curves  it  can  be  seen  that  wave  drag 
is  predicted  fairly  well  by  the  universal  wave  drag  function 
in  most  cases.   In  almost  every  case,  the  humps  and  hollows 
are  reflected  in  the  predicted  wave  drag  curve.   The  earlier 
rise  in  the  wave  drag  curve  for  high  values  of  C   and  the 
later  but  steeper  rise  in  the  wave  drag  curve  for  low  values 
of  C  are  both  well  predicted  by  the  universal  wave  drag 
function. 

The  wave  drag  of  vessels  with  lines  differing  greatly 
from  the  Taylor  Standard  Series  is  not  predicted  as  well  as 
the  wave  drag  of  vessels  in  the  Taylor  Series.   Some  sample 
curves  of  predicted  wave  drag  and  wave  drag  from  model  tests 
are  provided  for  comparison  in  Appendix  D.   In  most  cases, 
wave  drag  is  overpredicted  for  values  of  V//L  above  1.0  and 
underpredicted  for  values  of  V//L  below  1.0.   There  are  several 
possible  reasons  for  these  errors. 

At  high  speeds,  vessels  with  wide  sections,  full  water- 
lines  and  flat  buttock  lines  aft  will  derive  more  dynamic  lift 
than  vessels  with  deep  sections,  hollow  waterlines  and  full 
buttock  lines  aft.   As  a  result,  the  former  will  tend  to  squat 
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more  than  the  latter  at  high  speed.   The  vessel  with  the  deep 
sections  and  full  buttocks  aft  may  have  less  residuary  drag  at 
low  speed  than  the  vessel  with  wide  sections  and  flat  buttocks 
but  will  have  more  drag  at  high  speed  due  to  the  additional 
wave  drag  caused  by  squat. 

The  Taylor  Standard  Series  parent  has  deep  V-shaped  sec- 
tions and  full  buttocks  aft.   The  survey  launch  and  the  5.5 
meter  yacht  "Antiope",  on  the  other  hand,  have  very  straight 
buttocks  and  wide  sections  aft,  providing  more  squatting  sur- 
face than  the  other  vessels.  Although  the  destroyer  and  patrol 
boat  models  have  V-shaped  sections  aft,  they  are  not  very  deep 
and  the  buttocks  are  not  badly  curved.   Prediction  formulae 
based  on  the  Taylor  Series  would  therefore  be  expected  to 
underpredict  the  resistance  of  the  latter  three  vessels  at  low 
speed  and  to  overpredict  at  high  speed.   It  may  be  seen  from 
Appendix  D  that  the  patrol  boat  is  the  only  one  of  these  four 
which  does  not  behave  completely  as  expected. 

In  addition  to  the  differences  in  wave  drag  coefficient 
due  to  dynamic  lift  from  the  afterbody  lines,  other  differences 
in  the  lines  may  also  cause  errors  in  the  prediction  of  wave 
drag. 

Variations  in  model  size  and  testing  techniques  may  also 
explain  some  of  the  error  in  predicting  the  wave  drag  coeffi- 
cient of  those  models  sampled  which  were  not  in  the  Taylor 
Series.   As  discussed  previously,  flow  separation  and  wave 
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interactions  may  cause  scale  effects  over  a  small  range  of 
model  size.   Reference  9  discusses  scale  effects  in  hydrody- 
namic  testing  in  great  detail.   Many  plots  are  presented  of 
predicted  resistance  coefficient  versus  speed  curves  to  pro- 
vide comparisons  between  larger  and  smaller  geosims.   It  was 
noted  that  even  for  widely  differing  hull  forms,  from  merchant 
ships  to  racing  yachts,  the  smaller  geosim  has  a  larger  re- 
sistance coefficient  at  low  speeds  and  the  larger  geosim  has 
a  larger  resistance  coefficient  at  higher  speeds.   The  cross 
over  point  is  a  little  above  a  V//L  of  1.0.   This  is  shown 
roughly  in  Figure  23.   If  the  wave  drag  coefficient  prediction 
curves  in  Appendix  D  are  compared  to  Figure  23,  it  will  be 
seen  that  regardless  of  the  causes,  the  errors  show  strong 
similarities  to  model  scale  effects.   The  predicted  wave  drag 
coefficient  curves,  based  on  the  regression  analysis  of  the 
Taylor  Standard  Series,  follow  the  general  trend  expected  for 
the  larger  geosim.   The  model  test  resistance  curves  for  most 
of  the  models  having  hull  forms  unlike  the  Taylor  parent 
follow  the  general  trend  expected  for  the  small  geosim,  even 
though  they  were  actually  12.5  to  16  feet  in  length. 

Since  the  change  in  the  drag  curve  predicted  from  the 
Taylor  Series  caused  by  afterbody  lines  which  derive  more  dy- 
namic lift  at  high  speed  is  somewhat  similar  to  the  change 
caused  by  testing  a  smaller  model,  i.e.,  higher  drag  at  low 
speed  and  lower  drag  at  high  speed,  then  it  would  be  very 
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difficult  to  tell  which  effect,  if  any,  was  contributing  to 
the  underprediction  at  low  speed  and  overprediction  at  high 
speed  of  any  smaller  models  which  have  flat  buttock  lines.  The 
survey  launch  may  be  in  this  category  although  the  model  was 
16  feet  long,  and  it  is  more  likely  that  the  flat  buttocks 
caused  the  wave  drag  to  be  overpredicted  at  high  speeds. 

Since  "Antiope"  was  24.3  feet  in  waterline  length  and  the 
destroyer  was  a  20  foot  model,  it  can  only  be  their  relatively 
straight  buttock  lines  which  cause  their  wave  drag  coefficients 
to  be  underpredicted  at  low  speed  and  overpredicted  at  high 
speed. 

Although  it  is  difficult  to  tell  from  their  lines,  the 
trawlers  may  derive  slightly  more  dynamic  lift  at  high  speed 
than  the  Taylor  models  which  could  explain  the  underprediction 
of  their  wave  drag  at  low  speed  and  overprediction  at  high 
speed.   In  addition,  the  trawler  models  were  12.5  feet  in 
length  and  may  suffer  some  scale  effect  when  compared  to  the 
20  foot  Taylor  models. 

Lines  drawings  are  provided  in  Figures  8  and  25  through  29 
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XII.   FORM  DRAG  CONCLUSIONS 


Plots  of  form  drag  coefficient, 


l/2pV2LJ  =  CSCR(7L  "  '5)v%  (63) 


as  a  function  of  C   for  three  values  of  C   and  three  values  of 

P  v 

B/T  are  provided  with  data  points  in  Figures  15,  16  and  17  to 

show  how  the  data  was  faired.   Similarly,  plots  of  form  drag 

coefficient  as  a  function  of  C   for  three  values  of  C   and 

v  p 

three  values  of  B/T  are  provided  in  Figures  12,  13  and  14. 
From  these  graphs,  it  is  evident  that  form  drag  is  very  nearly 
proportional  to  both  C   and  C  .   Although  data  for  the  B/T  of 
3.0  is  plotted,  it  was  not  used  to  derive  the  equations  for  F 
and  F,  for  the  reasons  previously  discussed  in  Chapter  V.  Form 
drag  was  therefore  assumed  to  be  linear  in  B/T  between  the 
values  of  2.25  and  3.75. 

The  resulting  equations  relating  F   and  F,  to  B/T  and  C 
are  as  follows: 


F   =  .0027  +  (C   -  .5) [.021  +  .02033(|-  2.25)]     (64) 


F±   =    14.0  +  (|  -  2.25)[4.0  +  20.7(C   -  .5)]         (65) 


The  form  drag  coefficient  is  then 
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T  2 
Cf  =  (CQ  +  C1)  (±-)  (66) 


where  C   is  defined  by  Equation  (38)  and  C,  is  defined  by 


Equation  (36) ,  or 


Cf  =  (F   +  FlCv)  (£)  (67) 


With  each  of  the  plots  comparing  predicted  wave  drag  to 
model  test  results  is  listed  the  predicted  form  drag  and  that 
of  the  model  test.   Figure  19  also  contains  data  points  which 
represent  the  value  of  F,  calculated  from  the  equation  for  F 
and  the  residuary  resistance  coefficient  at  V//L  =  .5  from 
reference  7.   The  plotting  of  these  calculated  "data"  points 
for  F,  is  not  a  complete  demonstration  of  the  accuracy  of  the 
equation  since  it  assumes  the  accuracy  of  the  equation  for  F  . 
The  best  evaluation  of  the  accuracy  of  the  equations  for  F 
and  F,  is  probably  the  direct  comparison  of  the  predicted  value 
to  the  value  of  the  residuary  resistance  coefficient  at  a  V//L 
of  .5  from  reference  7. 

Form  drag  is  predicted  reasonably  well  for  all  vessels  in 
the  Taylor  Standard  Series  as  may  be  seen  from  Appendix  C. 
Prediction  is  best  for  those  vessels  with  volumetric  coeffi- 
cients in  the  middle  of  the  range  of  data,  approximately  .004. 

For  heavier  vessels,  form  drag  is  underpredicted  increas- 
ingly with  volumetric  coefficient.   This  effect  is  simply  the 
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manifestation  of  the  slight  nonlinearity  of  form  drag  with 

displacement.   The  effects  of  C   and  B/T  seem  to  be  adequately 

represented  by  the  Equations  (64)  ,  (65)  and  (66)  for  the  ships 

of  medium  and  heavy  volumetric  coefficients. 

For  extremely  fine  vessels,  the  effects  of  C   become 

P 

slightly  overpredicted.   As  a  result,  form  drag  for  vessels 

with  low  C   and  low  C   is  underpredicted  but  is  overpredicted 

for  vessels  with  high  C   and  the  same  C  .   This  is  another 

p  v 

manifestation  of  the  slight  nonlinearity  of  form  drag  with 
displacement. 

For  hull  forms  differing  from  the  Taylor  lines,  presented 
in  Appendix  D,  it  appears  that  form  drag  is  predicted  relatively 
well  for  slender  hulls  but  not  so  well  for  vessels  with  higher 
volumetric  coefficient.   The  old  destroyer  in  Appendix  D,  for 
example,  is  predicted  reasonably  well  with  a  volumetric  coeffi- 
cient of  .00139.   The  round  bottom  patrol  boat  in  Appendix  D  is 
predicted  less  accurately  with  a  volumetric  coefficient  of 
.00246.   The  survey  launch  in  Appendix  D  is  underpredicted  by 
about  500%  with  a  volumetric  coefficient  of  .0043.   Both  trawler 
models  are  underpredicted  by  over  100%  with  volumetric  coeffi- 
cients of  .00691  and  .00830.   The  form  drag  coefficient  of  the 
5.5  meter  yacht  "Antiope"  with  a  volumetric  coefficient  of 
.00627  is  predicted  reasonably  well  as  long  as  the  fairbody 
beam  to  draft  ratio  is  used.   The  form  drag  coefficient  predicted 
by  the  computer  program  for  "Antiope"  was  inaccurate  because  the 
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jfull  draft  was  automatically  used  to  calculate  the  beam  to 
draft  ratio.   A  misleading  beam  to  draft  ratio  of  only  1.35 

was  then  entered  in  the  prediction  formulae  for  F,  and  F  .   It 

r  1      o 

is  interesting  that  "Antiope"  is  predicted  so  well  when  the 
other  vessels  of  similar  volumetric  coefficient  are  so  badly 
underpredicted . 

Taylor  Series  models  with  high  volumetric  coefficient  are 
also  predicted  much  more  accurately  than  the  trawlers.   Few 
commercial  vessels  whose  feasibility  depends  on  carrying  capa- 
city and  acquisition  cost  can  justify  the  fine  lines  of  the 
Taylor  Standard  Series.   Merchant  ships,  for  example,  almost 
always  have  some  parallel  middle  body  and  shoulders  which  may 
generate  additional  wave  systems  and  induce  flow  separation. 
The  trawlers  have  propeller  apertures  and  may  not  have  after- 
body lines  as  smooth  as  the  Taylor  models  or  the  5.5  meter 
yacht.   The  Taylor  parent  hull  form  is  mathematical  and  care- 
fully avoids  harsh  curvatures.   It  is  difficult  to  conceive  of 
any  commercial  hull  form  which  would  treat  the  boundary  layer 
more  gently  than  one  with  the  Taylor  afterbody  lines.   A  5.5 
meter  yacht  is  a  racing  machine,  the  design  of  which  is  not 
constrained  by  practical  considerations  such  as  propeller  aper- 
tures and  carrying  capacity.   It  might  be  expected  therefore 
that  a  racing  yacht  would  have  less  form  drag  than  a  fishing 
trawler  of  similar  volumetric  coefficient.   Lines  drawings  of 
these  vessels  are  provided  for  comparison  in  Figures  25  through 
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29.   The  heavier  of  the  trawlers,  it  should  be  pointed  out,  is 
completely  outside  the  range  of  volumetric  coefficient  sampled 
by  the  regression  analyses  in  this  thesis.   If  the  lines  of 
the  trawlers  do  not  treat  the  boundary  layer  as  gently  as  the 
5.5  meter  yacht  or  the  Taylor  Series  models,  then  it  would  take 
only  a  small  region  of  separated  flow  to  explain  the  100%  dif- 
ference in  form  drag  coefficient  between  the  trawlers  and  the 
other  vessels. 

The  500%  error  in  predicting  the  form  drag  of  the  survey 
launch  cannot  be  explained  by  subtle  changes  in  the  afterbody 
lines  or  by  any  inaccuracy  in  the  prediction  formula.  With  a 
volumetric  coefficient  of  .0043,  a  reasonably  slender  vessel 

well  within  the  range  covered  by  the  prediction  formula,  a 

_3 
form  drag  coefficient  of  2.5  x  10    is  about  five  times  too 

large  and  indicates  that  the  boundary  layer  must  have  separated 
from  the  afterbody  or  else  the  transom  must  have  been  submerged 
deeper  than  the  design  waterline. 

Although  the  equations  for  F   and  F,  are  only  artificial 
formulas  fitted  to  the  data  rather  than  based  on  physical  laws, 
there  are  some  interesting  characteristics  to  the  gross  behav- 
ior of  F   and  F, . 
o      1 


F   =  .0027  +  (C   -  .5) [.021  +  .02033(|  -  2.25)] 


F1   =  14.0  +  (|  "  2.25)[4.0  +  20.7(C   -  .5)] 
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For  ships  with  low  values  of  C   near  .5,  skin  friction  error 
becomes  a  constant  independent  of  beam  to  draft  ratio.   For 
ships  with  low  values  of  beam  to  draft  ratio  near  2.25,  on  the 
other  hand,  F,  becomes  a  constant  independent  of  prismatic  co- 
efficient. 
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XIII.   WAVE  DRAG  RECOMMENDATIONS 

The  range  of  data  used  in  this  thesis  is  limited  to  that 
which  was  available  in  the  time  available.  The  1050  ship-speed 
combinations  generated,  while  a  substantial  number,  included 
speed  length  ratios  only  through  1.2.   It  is  recommended  that 
this  effort  be  extended  through  higher  speeds.   If  higher 
speeds  are  attempted  using  the  Taylor  Standard  Series  as  a 
data  base,  it  will  be  found  that  the  range  of  C   for  which 
data  is  available  is  narrowed  considerably  from  that  of  the 
lower  speeds,  eliminating  the  fuller  ships.   It  is  likely, 
therefore,  that  any  investigation  of  higher  ship  speeds  will 
have  to  be  done  using  either  destroyer  type  hull  forms  such  as 
Series  64  or  a  high  speed  Trawler  Series. 

In  addition  to  the  limits  on  speed  length  ratio,  the  1050 
ship-speed  combinations  generated  for  this  thesis  included, 
for  most  speeds,  only  those  values  of  C   between  .001  and  .006 
This  range  of  C   is  equivalent  to  a  range  of  displacement 
length  ratio  of  28.6  through  171.4,  i.e.,  very  slender  ships. 
These  values  of  displacement  length  ratio  are  limited,  for  the 
most  part,  to  military  vessels.   Extending  the  range  of  dis- 
placement length  ratio  to  400  or  C   to  a  value  of  .014  would 
include  most  interests  in  the  field  of  naval  architecture. 

If  the  sectional  area  curve  is  normalized  to  1.0,  then 

the  displacement  is  non-dimensionalized  and  becomes  equal  to 

C  .   If  L  =  1.0  and  A  =1.0,  then: 
p  x 


x=L 

f 


V  = 


A(x)dx  = 
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K)-dx  =  C   =  y  (68) 


x=0 


A  L  —     p    -o 


where  y   is  the  zeroth  moment  of  the  sectional  area  curve. 
The  longitudinal  prismatic  coefficient  may  be  thought  of  as  a 
zeroth  moment  of  the  sectional  area  curve.   Similarly,  a  first 
moment  of  the  sectional  area  curve  may  be  defined  as: 


x=L         x=l 

*  /  \  ^     f   x  A  ( x ) 
xA(x)dx  = 


dx  =  (C  ) (LCB)  =  y   =  y  LCB      (69) 


A  L  ~"    v  p'  %™->    Hl    poJ 
x=0  x=0    x 

and  LCB  becomes  the  ratio  of  the  first  moment  to  the  zeroth 
moment  of  the  sectional  area  curve.   The  second  moment  will 
then  describe  the  fullness  of  the  ends  more  specifically.   The 
third  moment  will  describe  the  asymmetry  of  the  sectional  area 
curve  in  greater  detail.   The  second  and  third  moments  together 
may  be  used  to  describe  the  momental  skewness  of  the  sectional 
area  curve.   These  and  higher  moments  of  the  sectional  area 
curve  may  be  defined  and  used  as  higher  order  terms  in  regres- 
sion analyses.   Weighted  variables  could  be  defined, as  was 
done  in  this  thesis, to  eliminate  the  dependence  on  beam  to 
draft  ratio.   With  a  sufficient  number  of  terms,  it  would  seem 
that  the  hull  form  could  be  adequately  described  so  that  wave 
drag  could  be  predicted  even  for  vessels  not  geometrically 
similar. 

In  this  thesis,  only  the  gross  first  order  effects  of 
displacement  and  beam  to  draft  ratio  were  collapsed  out  of  the 
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data.   Many  higher  order  effects  such  as  LCB  were,  by  the  use 
of  the  Taylor  Series,  carefully  avoided.   Although  the  wave 
drag  of  a  hull  should  be  predictable  from  a  few  well  chosen 
gross  parameters  without  addressing  every  detail  of  the  hull, 
there  are  many  variables  evaluated  in  regression  analyses  of 
total  drag  [16  and  23]  which  change  drastically  in  relative 
importance.   Reference  1  discusses  the  relative  contributions 
of  these  parameters  in  great  detail  and  presents  them  all  in 
a  single  compact  illustration,  re-drawn  as  Figure  24.   From 
the  behavior  of  the  correlation  of  these  parameters,  it  may  be 
anticipated  that  the  variables  will  have  to  be  very  carefully 
chosen  in  order  to  take  into  account  differences  in  hull  lines 
which  may  have  the  same  basic  form  coefficients. 
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NOTATION  FOR  FIGURE  2  4 


A,  =  total  section  area  at  the  AP 
A 

A  =  total  section  area  at  the  FP 

AP  =  after  perpendicular 

A  =  total  section  area  at  the  station  of  maximum  area 
x 

B  =  beam  at  the  station  of  maximum  section  area  on  LT,T 
x  WL 

CP  ■  v/(lwlV 

Sit  rtlwl/av2 

C   =  A  /(B  T  j 
X      X'    X  X 

f A  ■  VAX 

fB  =  VAx  : 

FB  =  position  of  the  longitudinal  center  of  buoyancy  measured 
from  the  FP 

FP  =  forward  perpendicular 

iR  =  angle  of  the  buttock  at  1/4  T„  at  station  L^  /20  forward 
of  the  AP  Wij 

i   =  entrance  angle  measured  between  the  centerline  and  a 
tangent  to  the  designer's  waterline  at  the  FP 

L  =  ship  or  model  length  measured  on  the  designer's  waterline 

R™  =  total  resistance  of  the  hull  in  pounds 

S  =  wetted  surface 

(S)  =  S/V2/3 

T  =  draft  measured  from  the  designer's  waterline  to  the  lowest 
point  of  the  station  of  maximum  area 

V  =  ship  speed  in  knots 
A   =  displacement  in  tons 

V  =  volume  of  displacement  in  cubic  feet 

4  =  A/  ( .  01LTTT  )  3 ,  displacement-length  ratio 

Till 
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31.56  ft 

T.W.L. 

24.3  ft 

Beam 

6.27  ft 

Draft 

4.64  ft 
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XIV.   FORM  DRAG  RECOMMENDATIONS 

Due  to  time  limitations  and  to  the  relatively  small 
contribution  of  form  drag  to  the  residuary  drag  at  high  speeds, 
a  computer  regression  analysis  was  not  utilized  to  derive  for- 
mulae for  F   and  F, .   The  functional  forms  of  F   and  F,  were 
o      I  o      1 

instead  determined  graphically  and  were  limited  to  linear 
functions  of  C   and  B/T.   Although  these  simple  linear  predic- 

ir 

tion  formulae  work  surprisingly  well,  there  are  some  data 
points  plotted  after  the  fact  which  are  not  well  predicted.  A 
regression  analysis  sampling  all  315  residuary  resistance  co- 
efficients at  a  V//L  of  . 5  or  less  from  reference  7  would 
surely  pull  in  most  of  the  wild  points,  particularly  if  higher 
order  terms  were  made  available  for  stepwise  inclusion. 

If  a  regression  analysis  were  set  up  to  fit  prediction 
formulae  to  form  drag,  it  would  take  only  a  little  additional 
programming  to  include  other  higher  order  effects  such  as  LCB 
location  in  the  analysis. 
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APPENDIX  A 


Program  Used  for  Calculation  of  Weighted 
Variables  and  F~  for  the  1050  Ship-Speed  Combinations 
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1 
~2 

3 

4 

5 

— 6~ 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
-17 


AN  ANALYSIS  OF  RESICUARY  DRAG 


DIMENSION  XSLGP( 11.11) ,ZSLCP( 11,11 ),OAl 11,11) 

DIMENSION  YY(  ll,ll),Y(ll,il)  ,AACH,li),A(llJ  ,/JM{ll),G(ll) 

DIMENSION  CR(15f7) 

INTEGER  VL,CVV,8C 

REAL  M 

REAL"  LCB 

REAL  LBP 
REAL  LWL 
INTEGER  WL 
INTEGER  STA 
BUG*. 5 

LWL=190.   "  

LBP=190. 

T=4.656 

8=17.461 


RE AD (5, I )     ( ( YY(WL.STA),WL=1,9) ,STA=1, 11) 
t-FORMAT  (  9F8  .  7  ) — 


18 
19 

20 


21 
22 
23 
24 
25 
26 
27 


RE AD (5, 39)     ( t CR ( VL ,CVV ) , VL= 1 , I  5 ) ,CVV= 1, 7) 
39    FORMAT    ( 15F5.2) 

DO    35    STA=l,ll 


ADDITION    OF    WATERLlNES    .9    AND    .7    TO    THE    TAYLOR    CFFSETS 
—    DO    20    1=1,7 

WL=12-I 

YYUL,STA)=YY(WL-2,STAI 

20    CONTINUE 

YY( 3,STA)=\Y(2,STA)    

YY(2,STA)=(YYd,STA)*YY(3,STA)  )*0.5 

YY( 4,STA)=(YY(3,STA)+YY(5,STA) )*0.5  


28 

29 

-30 


31 
32 
33 
34 


35 

36 
37 
38 
39 
40 
41 


35    CONTINUE  .  

CVV=1 
t5C0N  T I NUE 

PROPORTIONING    OF    FULL    OFFSETS    FRCM    TAYLOR    FRACT  ICNAL    OFFSETS 
DO    7    STA=l,ll 
DO    7    WL=l,ll 

Y(WL,STA)=YY(WL,STA)*B*0.5 
—  7   CONTINUE  " 

CALCULATION  OF  SLOPES  AT  EACH  CFFSET  REFERENCE  PCINT 
DO  21  WL=1,11 
DO  21  STA=2,  10 

21  XSL0P!WL,STA)=<Y(WL,STA*l)-YlWl,STA-l))*05.0/LWL 

DO  22   STA=l,il 

DO  22    WL=2,1C 

22  ZSLOP(WL,SrA)=(YIWL+l»SrA)-Y(HL-l,STA))*05.0/T 
DO  23  WL=l,ll 
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2  23   XSLOP(WL,l)=2.0*XSLCP(WL,2l-)<SLCP(WL,3) 

3  DO    24    WL=l,ll 

4  24    XSLCPCWL,  U)=2.0*XSL0PIWL,  IC l-X SLOP ( WL , 9) 

-5 00    25    STA=I,ll 

6  25    ZSLCP( l,STA)=2.0*ZSLOP(2tSTA)-2«LQP(3fSTA) 

7  00    26    STA=1,  11 

>8  26    ZSLOPI  U,STA)=2.0*ZSLCP(10,STM-ZSLCP<9,STA) 

9  VL»1 


il  BUGR=BUG*-5 

c 

>2  WRITE    (6,191     flUGfi 

»3  19    FORMAT     (•  l» , •♦**♦••**»•*»••*•••»*•*••*****•       SHIP    NUMBER    »     «,F5.2, 
C '        *«*♦*****♦»♦»»****»*»*»♦*•*♦#*••♦*»«**•»♦*»*»»*♦*♦*»*•) 

»4  FR-0.0 

>5  GO    TO    34 

C 

>6  6    V*SL*SQRTILWL)*1.689 

1 C CALCULATION  OF  EXPONENT  FOR  WEIGHTING  INCREMENTAL  AREAS  " 

>7  FR»32. 2/(V**2) 

C 

C  CALCULATION    CF    FRQUCE    NUMBER 

58  F«V/U32.2»LWL)»«.5)  " 
C 

59  34  CONTINUE        " ~ — 

SO  S  T  A  =  I 

61  2  CONTINUE 
C 

C  INTEGRATION  OF  WEIGHTED  INCREMENTAL  AREAS  OVER  THE  DRAFT  FOR  EACH  SECTION 

62  DO  31  wL=l,ll 

63  31  AAlWLtSTA)  *YlWL,STAJ*EXP(-FR»IH-ll»T/lO.O)*2.0 

64  EVNN=0.0 

65  OOOC*0.0 

66  00  32  Wl  =2,  10,2 

67  32  EVNN=EVNN*AA(WL, STAJ 

68  DO  33  WL*3,9,2 

69 33  OOOD=OODO*AA(WL,STA)  " — 

70  A(STA)=(T/30.0)*(AA(1 , ST A ) *4 ,0»E VNN ♦2.0*0C0D ♦ I A ( 1 1 , ST  A ) ) 
C 

C  CALCULATION  OF  LONGITUDINAL  MCfENTS  OF  SECTIONAL  AREAS 

71  AMISTAI=AI STAI*FLCAT(STA) 
C 

72 — — IFIFR.GT.0.0)  GO  TO  18  — 

C 

C  INTEGRATION  OF  INCREMENTAL  GIRTHS  OVER  EACH  HALF  SECTION 

73  DO  27  WL=ltll 

74  2  7    DAI  hL,STA)=SQKT(  1  .0*1  XSLOPC  WL  t  ST  A)  I  **2*  I  ZSLCP  IWL  t  ST  A}  )**2) 

75  EVN-0.0  

76 ODOC'O.O  — 

77  DO    28    WL=2,10,2 

78  28    EVN*EVN*DAIWL,STA) 

79  DO    29    WL»3»9,2 
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80  29    ODDOODOO*DA(HL,STAI 

81  G{STA)=(T/30.0)*IDA(l,STA)»4.0*EVN*2.0*G00D»0*(lltSTA) 1 

C 

82 1/8  TONTINUE 

C 

83  STA*STA*l 

84  IF     (STA.LE.U)    GC    TC    2 
C 

C  INTEGRATION    OF    SECTICNAL    AREAS    OVER    LENGTH    FOR    DISPLACEMENT 

85 000  =  0.0  — • • 

86  EvEN=0.J 

87  00  3  STA=2t 10,2 

88  3  EVEN=EVEN*A( STA) 

89  DO  A  STA=3,9,2 

90  4  000=000«-A(  STA) 


--9 1 OlSPL=(LBP/30.IMAm*<».*EVEK#2v*0DD*Ariiri 

C 

C    INTEGRATION  OF  MOMENTS  OF  AREAS  CVER  LENGTH  FOR  LCB 

92  AT=C.O 

93  AHT=0.0 

94  DO  9  STA  =  1,H 

"95   9  AT»AT*A(STA)        — —  — — — — 

96  .         DO  5  STA=l,ll 

97  5  AMT*AMT*AM(STAl 

98  LCB*lAMT/AT-6.)*ILBP/LWl)*.l 

99  M»LhL/<DISPL**. 333333333) 

C 

100  CP«0ISPL/CAC6)*LWL) 

101  CPS=CP 
C 

102  CV=»0ISPL/(LWL**3» 

103 CVST=CV*1000.0 " 

C 

104  CM*A(6)/(B*T) 
C 

105  CB=CM*CP 
C 

^06 IF    tFR.LE.  O.VZi    GCTTOl* ' — ~ — 

C 

C         F2    FOR    CW    BASED    CN    L**2 

107  CCR=CR(VL,CVV) 

108  CR5=CR( I ,CVV) 

109  CR(VL,CVV)=    0.001*CRl VL  ,CVV) 

110 F2»ICH(VL,CVV»  "-    CRIl.CVV))     ♦   ^L2   *    C(l LWL**3 1/ (OISPL ) » **2> 

C 

111  PUNCH    «<►     ,     SL.CCP,CR5,CVST,CVT  ,CPS,CPN,SL2,BOR,F2 

112  44    F0RMATCF«.2,2X,F4.2.2X,F4.2,2X,F5.3,2X,F5.3,2X,F6.<»,2X,F6.4,2X, 

CF6.4,2X,F<>.2,2X,F6.2) 
C 

113 —WRITE     (6,13)     SL,F, OISPL, CP,CK,CB,CV,CCR, LCB, F2      — 

114  13    FORMAT  (  'O'  ,  F  4.  2  »  2  X  ,F5  .  3  .  3X  ,  F  fc  .  C  .  3X  ,  F9.<»  ,  3  X  ,F9  .<.  ,3X  ,  F9.4  ,3X  ,F9.6  ,  3X 

C,F6.«,3X,F9.4,5X,F12.6) 
C 
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115  SL»Sl>.05 

116  VL»VL*1 

117  IF  (SL.LE. 1.2)  GO  TO  6 

— 118 GOT  TO  36 

C 

119  14  CONTINUE 
C 

C  INTEGRATION    OF    GIRTHS    OVER    LENGTH    FOR    WETTED    SUPFACE 

120  EVE=0.0 

_12l     OD*  CO  • — — 

122  00     16    STA=2.10,2 

123  16    EVE=EVE«-G(STA) 

124  00     17    STA=3,9,2 

125  17    OD=CD*Gl STA) 

126  S=0.06667*LBP*(G(l)*4.«fcVE*2-*CC*G(  ID) 


C         CORRECTION    TO    WETTEC    SURFACE    FCR    SKEG 

127  S=S-T*LWL*(0. 3372-0. A38*CP) 
C 

128  WS*S/SQRTtDISPL»LWL) 
C 

129     CPN=CP " 

130  DISP=DISPL 

131  CVT=CV*1000.0 

132  VR=CV 

133  SL2  =  S/(LWL»*21 

134  BDR*B/T 

135 F0=  .0027MCP-.5)  M .021 ♦ I 30R-2. 25)*. C203331 

136  F0=F0*0.001 

137  Fl*  !«.♦ (6DR-2.25)*(4.»20.7*(CP-.5) ) 

138  Fl*Fl*0.001 

139  FF«VR*Fl    ♦    FO 
C 

140  ~        WRITE(6, ll)LWL»B,T,OISPL,CP»CH  ~ 

141  11    FORMAT! •-• ,»LwL     ■• t F7-2 t 5X» • 8    =     •,F6.2,5X,'T    =     • , F6.2 ♦ 5X, • 01 SPL    = 

C     •,F8.l,5X,,CP    =     •,F9.5,5X,»C>'    =    «,F9.5) 
C 

142  WRITE     (6,12)     CB.CV.S.WS 

143  12  FORMAT! •-• ,'CB  =  •  ,F S. 5 , 5X , • CV  =  ',F9.5,5X,'S  *  • ,F 10. 0, 5X , • CS  -  • 
CiF9.4)      ; -<— 

C 

144  WRITE(6,37)  SL2tFF 

145  37    FORMAT  I'-'  ,' S/LWL**2    =     • , F 1 0.8 , 10X , • FO    ♦     F 1 «G ISPL/ ( LWL **3)     ='tFl0. 

C81 
C 
1*6     WRITE  (6,  8) 

147  8    FORMAT! •-•, 'SLR'  ,5X,'F«  ,5X, 'CISFL     *«,5X,,CP    ♦,,8X,'CM    *«,9X,,CB    *• 

Ct9X,»CV    *,,5X,»       CCR  '.SX.'LCE    ••,10X,•F2,) 

C 

148  GO    TO    6 

C  

-149      "36   CONTINUE  

150  CVV=CVV«-l 

C  DIMENSIONS    OF     A00IT1CNAL     SHIPS    TC    BE    CALCULATED 

151  BUG*BUG»1.0 
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152  T=6.5851 

153  B=24.694l 

154  IF  (BUG. LE. 2.0)  GG  TO  15 

155  T=8.065 

156  6=30.244 

157  IF  (BUG. LE. 3.0)  GC  TO  15 

158  T=9.3127 

159  B=34.9227 

160  IF  (BUG. LE. 4.0)  GG  TO  15 

161  T=10.4ll9 

162  B=39.0447 

163  IF  (BUG. LE. 5.0)  GC  TO  15 

164  T=ii.406 

165  B=42.77 

166  IF  (BUG. LE. 6.0)  GC  TO  15 

167  T=12.3196 

168  6=46.1984 

169  IF  (BUG. LE. 7.0)  GC  TO  15 
C 

170  WRITE(6,100) 

171  100  FORMAT! • V t • END  CF  PROGRAM') 
C 

172  STOP 

173  END 

SENTRY 
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APPENDIX  B.l 

Listing  of  Control  Cards  and  Some  Output  Examples 
for  One  Regression  Run 
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Sxample  of  first  page   of  printout   for  V/7L  =  1.0  run 


rATISTICAL    PACKAGE    FOR    TH  F    SOCIAL    SCIENCES    SPSiiH    -    RPLRASE    6.02 


03/25/76 


SPACE  ALLOCATION  FOR  THIS  RUN.. 
TOTAL  AMOUNT  REQUESTED 


80000  BYTES 


DEPAULT  TRANSPACE  ALLOCATION 


10000    BYTES 


MAX     NO    OP    TRANSFORATIONS    PERMITTED 

HAX     NO    OF    RECODE    VALUES 

MAX    NO    OF    ARITHH.OR    LOG . OPERATIONS 


fiTSULfl  NG    WORK  SP ACE~SLtOCAf ION 


100 
U00 
800 


"TOO  00    BYTES 


RUN    NAME 
FILE    NAME 
VARIABLE    LIST 
VAR    LABPLS 


MULTIPLE    REGRESSION     ANALYSIS    O?    THE    TAYLOR    SERIES 
DATA     FOR     PREDICTION     OP    WAVE    DRAG 
SLR,CR,CR5,CVST,CVT,CPS1,CP,SL2.BDR,P2 
SLB     _SPFFC    LENGTH     3ATIO    / 

TAYLOR    RESIDUARY    DRAG    COEFFICIENT    TIN  E5~TiJ0U    / 
CR    AT    A   SLR    OP    0.5      / 

WEIGHTED    VOLUMETRIC    COEPPICIENT    TIMES     1000   / 


Cfl 

CR5 

CVST 

CVT 

CPS1 

CP 

512 


ACTUAL    VOLUMETRIC    COEFFICIENT    TIMES     1000    / 
WEIGHTED    PRISMATIC    POLYNOMIAL    TERM    OF    ORDER    ONE    / 

ACTUAL    PRISMATIC    COEFFICIENT    /  

"TfX'TrO   OP    WETTEDSURPACE    TO    LENGTH   SQUARED"-/" 


INPUT    FORMAT 


BDR  BEAM    TO    DRAFT    RATIO 

P2  CV    STAR    SQUAR2D    TIMES    CW    BASED    ON    L       / 

FIXED    (PH.2,2X,F4.2,2X,?!».2,2X,P5.3,2XrF5.3,2X,P6.U,2x. 

P6.U,2X,F6.4,2X,P4.2,2X,F6.2) 


ACCORDING  TO  YOUR  INPUT  FORMAT,  VARIABLES  ARE  TO  BE  BEAD  AS  POLLOWS 
BECORD 


VARIABLE 

FORMAT 

SLR 

p  a.   2 

CR 

F   U.    2 

CR5 

p  a.   2 

cvsr 

P    5.    3 

CVT 

P    5.    3 

CPS1 

P    6.    4 

CP 

p   6.  a 

SL2 

T   6.    if 

BDR 

p  y.  2 

P2 

P    6.     2 

HE  INPUT  PORMAT  PROVIDES  POR   10  VARIABLE 
T  PROVIDES  FOR    1  RiCOHDS  ('CARDS')   PER  C 


COMPUTP 

COMPUTE 

COMPUTE 

COMPUTE 

COMPUTP 

COMPUTE 

COMPUTE 

N  OP  CASES 

REGRESSION 


STATISTICS 


COLUMNS 

1- 

0 

7- 

10 

13- 

16 

19- 

23 

26- 

30 

33- 

38 

m- 

U6 

T»9- 

"5U 

57- 

60 

63- 

68 

SE. 


10    WILL    B"5    READ 
A     MAXIMUM    OP         68    'COLUMNS'     ARE    USED    DN    A    RECORD. 


P3»P2*  (CVST** (0.2)) 

CPSV*  (CPS1-.UR)/.?7 

CPSA-SIN  (3.  1U*CPSV) 

CPSR=SIN(3.14*CP-;V«2.) 

CPSC  =  SIN  (3.  14*CPSV*  3.  ) 

CPSD=SIN  (3.  1<**CPSV*'I.) 

CPSE=SIN  (3.  1<4*CPSV*5.) 

66 

VARIABLES    =  P3,CPS1, 

CPSA,CPSB, CPSC,CPSD,CPSE       / 

REGRESSION-FJ    WITH    CP31(1), 

CPSA(  1)  ,  CPSB(1)  ,CPSC(  1)  ,  CPSD  (1)  ,CPSE(1)  t 

RESID=0    / 

ALL 
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Example   first  page   of  printout   for  V//I  ■  1.2  run 

STATISTIC*!    PACKAGE    FCR    THE    SOCIAL    SCIENCES    EESSH    -    RELEASE    6.02 

0  V2S/76 

SPACE    ALLOCATICN    FOR    THIS    BUN.. 

TOTAL    AMOUNT    REQUESTED                                                                           80000    BYTES 

........ 

DEFAULT    TPANSPACE     ALLOCATION                                                           10000    ATTPS 

MAX    NC    OP    TRANSFORMATIONS    PERMITTEE                     100 
MAX    NO    OP    RF.CODE    VALUES                                                      «00 
MAX    NO    OP    ARITHM.OR    LOG .OPERATIONS                       800 



T?5Tn;TTNG~TOBirSTAXT~AtrOCTTTC> TOXJaO^BTTP^ 

'RUN    NAME  MULTIPLE    REGRESSION    ANALYSIS   OF    THE    TATLOR    SERIES 

FILE    NAME  DATA    FOR    PREDICTION    OP    WAVE   DRAG 

VARIABLE    LIST  SLR ,CR ,CR e ,CV 5T, CV? ,C PS  1, CP, SL2, BDP, P2 

VAR    LABELS  SLR  SPEED    LENGTH    RATIO    / 


~Cfi  TAYlOR    RESIDUARY  ~DRAG"CT£FFICrENTTrNES    1000    / 

CR5      C9    AT    A    SLR    OP    0.  S      / 

CVST         WEIGHTED    VOLUMETRIC    COEFFICIENT    TI"ES    1000     / 
CVT    _    ACTUAL    VOLUMETRIC    COEFFICIENT    TTMFS     1000    / 
fflB'1     '    WEIGHTED    PRISMATIC    POLYNOMIAL    TERM    CF    ORDER    ONE 
CP         ACTUAL    PRISMATIC    COEFFICIENT    / 


ST2         RATIO- OP    WETTED"~SURFACE    TO    LE  S'5Tff  "  STC0  A  RETT"/ 

BDR         EEAM    TO    DRAPT    RATIO 

P2  CV    STAR    SOUAREE    TI«PS    CW    EASFD    ON    L       / 

^NPUT   POPMVT FIXED     (P«  .  2  ,  2X  ,  ?U  .  2  ,  2X  ,  pa .  2,  2X,  F  5 .  3,  2X,  F5  .  3,  ?X ,  P6  . 1 ,  2X, 

F6.1»,2X,F6.a,  2X,F«.  2,2X,?6.  2) 

ACCORDING  TO  TOU8  INPUT  FORMAT.  VARIABLES  ARE  TO  BE  J*FAD  AS  FOLLOWS^ 
VIBIAELE   FORMAT   BECCRC     COLUMNS 


SIR  F    tt.    2  1  1- 


en nr;~2 1 r^—ro 

CB5  F    U .    2  1  13-16 

CVST  " f    5.     3  1  19-      23 

CVT  F   5.     3  1  26-       30 

CFS1  F"    6.    0  1  33-      38 

Cf  P   6.    U  1  11-       U6                                                                                                 

" "SIT F    6.    a 1 H  <n--  ~5H 

FDR  F    U.     2  1  57-       60 

F2  f   6.    2  1  63-      68 

THE    INPUT    FORRAT    PFOVIEFS    POR       10    VARIAELES.  10    ait  I    RE    BPAD 

IT    PR0VICE3    POR       1     RECORDS     ('CARDS')     PER    CASE.        A    MAXIMUM    OP         68     'COLUMNS'    ARE    USED    ON    A    RECORD. 

COMPUTE  F3  =  F2*  (CVST**  (0.2) )  

CCFPUTE  CPSV=(CS1-.M)/.27 

COMPUTE  CPS A=SIN (3. 1U*C?SV) 

COMPUTE  CPSD=STF(3. 1U*CPSV*2.) 

COMPUTE  CPSC=SIN (3. 1«*CPSV*3.)  

COMPUTE CPSD=J5IN(3.JU*CPSV*<».) 

COMPUTE  CPSF=SIN(3. 14*CPSV*57J 

N    CF    CAS?S  S8 

BEGRESSICN  VARIABLES    =P3,CPS1, 

CPSA,CP3B,CPSC,CPSD,CPSE  / 
REGRFSSTON=P3  WITH  CnS1(1), 
CPS  A  (1)  tCESB(1)  ,CPSC  (1)  ,CPSC(1)  ,CPSE(1)  , 


RESID=0    / 
STATISTICS  ALL 
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T?IB    SIGRFSSICS    ANALYSIS    OP    TRF    TAYLOR    SIBITS 

I     OATA (CREATION    CATE    »    C3/25/76)  FOR    PREDICTION    OF    WAV    DRAG 


0  3/25/76 


PAGE  10 


»    1    T    I    PL    E      REGRESSION 


F)ENT    VARIABLE:       P3 


TTBSE'fiTrir 
P3 

82. 19550 
67.CC98  1 
93.21245 


PFOI 


VAPIAELE    IIST 
REGRESSION    LIST 


83.2CC78 
81.48373 
82.57573 
8ft.  16908 
82.S1«0ft 
T5T99"«a7~ 
57.815H0 
58.93739 
62.5G374. 
68.24979 
68.35806 
"557301 
55.72318 
'Si. 26572 
57.18381 

~  58.92725 
62.78528 

~ 5^736  976- 
<.3.95481 

~~IT3. 73724 
U7.15891 

~  H9.1207C 
52.91765 

—isrirs'iis- 

40.51698 
"ft  1.9902 8 
44.05937 
45.86856 
48.  148e5 

""40778557' 
46.65015 

"  47.  17084 
48.31050 
48. 34575 
50.66  368 

71139a" 

.65862 
50.21512 
47.43613 
50.37614 
53.21028 
53.71008 
71.JP788 
73.35895 
7J.51 393 
74.56857 
77.351)0 
83.68127 
90.330  12 
e7. 89235 
8!. 33984 
83.1 1327 
85.CC676 
31.58952 


"ft?. 
52. 


"■p'RfiBic'mr 

F3 

84.96455 
85.21678 
85.41003 
~~ 85774049" 
86.02934 
eft. 17039 
84.5  897  1 
84.96455 
85.ft0820 
85. 85475" 
60.79886 

"60.39572 
59.93355 
59.41650 
58.78702 
357  OH 
62.25871 
61.69379 
61.07109 

"60.39572 
59.60905 

^5  8T60"T?3~ 
44.93106 

~IKT.78830~ 
44.62337 
ftft. 4565ft 
44 .22424 

"TJ  T9055TT- 
45. 38956 
ft5. 2120ft 
45.01193 
4ft  .78830 
44.51935 
44.  17999— 
48.49669 
48. 96924 
49.59085 
50.29916 
51.09856 

-5272777  0— 
46. 89206 
47.454  30 
48.  16522 
48.06924 
50.01  167 
51.46852 
81.*0811 
81.60899 
81.59581 
81  .51459 
81.505  7  1 
.39453 
,51022" 
.56429 
.60129 
, 60899 
81.57910 
81.47478 


81 

81, 

Rl 

81 

81, 


SPSIDUAL 

-2.769071 

1.7930  26 

7.7324C5 

~5778ft05B~ 

-2.828582 

-2.6666fft 

-2.0  1^988 

-.7954845 

-2.S9D177 

~T1392Tr75~ 

-2.983459 

-1.1158  3  32 

2.57C189 

8.8332P1 

9.571038 

;r«?6r 

•6.535548 
■5.428079 
■3.987284 
■1  .468479 
3.176223 


PLOT~OP    STASDASBTZ;D"~HT§lDnirir~ 
■2.0  -1.0  0.0  1.0  2.0 


~*r 


-.2  Uht  W 

-1.076247 

-1.07136  3 

2.535529 

4.664150 

9.6934P1 

-7.5C61ir 

•4.872562 

-3.221763 

-.9525656 

1.080257 

3.6295C3 

■7.753355" 

-1.846557 

-1.7984C4 

-1.280359 

-1.953417 

-.4348861 

-5.ft7626-r 

5.766555 

2.764902 

-.7291036 

2.007595 

3.1966C8 

-.7584471 

-8.2  202  51 

•B.?50Ce8 

-8.081997 

-6.996025 

-4. 1464'0 

2.2  867  18 

"9.11999ft 

6.329051 

3.738543 

1.50U259 

3.ft296tt0 

10.  11471 
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HOLTTPLE    REGRESSION     ANALYSIS    CP    THE    TATLCP    SEflES 

FILE  DATA (CREHIPH    CATF    =    03/25/7€) 

FOB    PREDICTION    OP    HAVE    DRAG 

*******       PLOT:       STANEARCIZEC    RESIDUAL     (EOH'J)     -• 

PREDICTED    STANDARDI7PD    DEPENDENT    VARIABLE     (ACROSS)        « 


EIPENDENI    VARIABLE:        F3 


VAfil  ABI  !    LI51        1 
REGRESSION     IIST        1 


-2.0                 -1.0                   0.0                    1.0 

2.0 

2.0    ♦                                                       I 

.    _ . —  — __ . 

i 

1.0    ♦                                                          T 

.  *.._ 

.           __  — _                     _  ._      .  —  .      .  _ 

.  —        j 

T                                                     *       I                            * 

I                               *                       *    I                            *    • 

I                                                             I                            * 

I                               **                          I                                 * 

—  -      i                       -  *   •*             2   I~               ~~"2 



I                              2       *                   I                        •*    ♦ 

I                              3    *:            *2    I                              * 

I                              2                      **I                           **3 

T                                 *          *               *  I 

I                                 *                           *I                              * 

I                                                                I                              * 

.                               j       __ _ 

T 

-1.0    ♦                                                             T 

-2.0    ♦                                                         T 

-2.0                -i.O                   0.0                   1.0 

2  .0 

BOWS,  COLUMNS    1:       VALUES    OUTSIDE     (-3.0,3.0) 

POWS,COLTJHNS    X:       VALtJES    TN     (-3.C-2.05)     OR     (2.05,3.0) 
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APPENDIX  B.2 

Listing  of  the  Resulting  Regression  Coefficients 
for  Each  of  14  Runs 


Coefficients  are  listed  in  column  "B1 
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F3 

a 

6.2C  30 

* 

a 

0.5400 

F3 

= 

3.8  82  5 

CP 

4c 

= 

0.C4C0 

F3 

= 

5.7565 

* 

= 

0.5450 

F3 

= 

3. 5908 

CP 

4c 

= 

0.c45  0 

F3 

= 

5.3567 

* 

= 

0.5500 

F3 

= 

3.  350  7 

CE 

4t 

= 

0.5500 

F3 

= 

5.C274 

* 

= 

0.5550 

F3 

= 

3. 1836 

CP 

* 

= 

C.c550 

F3 

= 

4.7863 

* 

= 

0.5600 

F3 

= 

3.  1060 

CF 

* 

= 

0.5600 

F3 

= 

4.6529 

4c 

a 

0.5650 

F3 

= 

3.1292 

CP 

4c 

= 

0.5650 

F3 

= 

4.6283 

* 

= 

0.5700 

F3 

= 

3.2584 

CP 

4c 

a 

0.5700 

F3 

= 

4 >  1 1  4  3 

4c 

a 

0.5750 

F3 

a 

3.4928 

CP 

4c 

= 

C.575C 

F3 

= 

4  .9041 

* 

= 

0.5800 

F3 

a 

3.8259 

CP 

4c 

= 

0.5800 

F3 

= 

5.  1854 

4c 

a 

0.5850 

F3 

a 

4.2463 

CP 

4c 

= 

0.5850 

F3 

= 

5.5415 

4c 

a 

0.5900 

F3 

= 

4 .7392 

CP 

4c 

a 

0.5900 

F3 

= 

5.9C3U 

* 

a 

0.5950 

F3 

= 

5.2877 

CP 

4c 

= 

0  .5950 

F3 

= 

6.40  16 

* 

= 
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= 

5.8743 

CP 

4c 

a 
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F3 

a 

6.P665 
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= 

0. 6050 
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6.4828 

CP 

4c 

a 
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F3 
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7.3398 

* 

= 

0.6100 

F3 

= 
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CP 

4c 

a 

C.6  1CC 

F3 

a 

7.8^*2 

4c 

a 
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F3 

= 

7.7134 

CP 

* 

a 
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F3 

= 

P. 2645 

4c 

= 
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F3 

= 

8.  3193 

CP 

4c 

= 
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F3 

= 

8.7177 

4c 

= 
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F3 

= 

8.9156 

CE 

4c 

= 
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F3 

a 

9.  1748 
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= 

0.6  300 

F3 

= 

9.5055 

CE 

4c 

a 

0  .6300 

F3 

= 

9.6495 

* 

= 

0.6350 

F3 

= 

10.0956 

CP 

4c 

= 

0.6350 

F3 

= 

1C.  1594 

4c 

= 

0.6400 

F3 

= 

10.6958 

CP 

4c 

- 

0  .6400 

F3 

= 

10.7230 

* 

= 

0.6450 

F3 

= 

11. 3170 

CP 

4c 

a 
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F3 

= 

1  1  .3579 

4c 

= 

0.6500 

F3 

= 

11. 9704 

CP 

4c 

= 

0.6500 

F  3 

= 

12.076^ 

4c 

= 

0.6550 

F3 

= 

12.6654 

CP 

* 

a 

^.6550 

F3 

= 

12.8928 

4c 

= 

0.6600 

F3 

- 

13.4088 

CP 

4c 

= 

0.66CC 

F  3 

s 

13.8021 

4c 

= 

0.6650 

F3 

- 

14. 2029 

CP 

4c 

= 

C.6650 

F3 

= 

14  .7979 

* 

= 

0.6700 

F3 

= 

15.0454 

CP 

4c 

= 

0.6700 

F3 

= 

15.8621 

* 

a 

0.6750 

F3 

= 

15.9286 

CP 

* 

= 

0.6750 

F3 

= 

16.9668 

4c 

= 

0.6800 

F3 

a 

16.8395 

CP 

♦ 

a 

0.68C0 

F3 

= 

IP. 0753 

4c 

= 

0. 6850 

F3 

= 

17.7606 

CP 

4c 

s 

C  .6850 

F3 

= 

19.  1438 

4c 

= 

0.6900 

F3 

= 

18.6708 

CP 

4c 

a 

0.6903 

F3 

a 

20.1239 

4= 

= 

0.6950 

P3 

= 

19.5468 

CP 

4c 

a 
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F3 

= 

20.96^3 

4" 

a 
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F3 

= 
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CP 

4c 

= 

C.7CCC 

F3 

= 
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CP    * 

= 

0.480  0 

F3 

- 

13.3153 

CP 

4c 

= 

0.4800 

F3 

X 

20.2021 

cp  * 

= 

0.4  85  0 

F3 

- 

13.6364 

CP 

4c 

■x 

0.4850 

F3 

= 

20.3639 

CP    * 

-= 

0.4900 

?3 

— 

1  3.90^6 

CP 
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3 

0.4900 

F3 

s 

20.4780 

CP    * 

— 

0.4950 

F  3 

- 

14.0911 

CP 

* 

B 

0.4950 

F3 

= 

20.5006 

CP    * 

= 

0.5000 

F3 

- 

1 4 ,  ]U U  u 

CP 
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F3 

X 

20.3951 
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= 

0.505  3 

F3 

— 

14.0436 

CP 

* 

= 
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20. 1353 
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D  .  5  1 0  0 
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1 3.77  5  2 

CP 

* 
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0.5100 

F3 
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= 
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= 

0.5300 

F3 

= 

1 1.2241 

CP 

* 

= 

0.5300 

F3 

= 

16.4848 

CP    * 

= 

0.5  3  50 

»3 

= 

10.  3668 
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F3 

= 
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F3 
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7.7  309 
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♦ 

3 
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F3 
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CP 

* 

x 
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F3 

= 
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cp 
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= 

11.4610 
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= 

0.6  20  0 

73 

- 
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CP 

* 
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F3 
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= 
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= 
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F3 
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= 

14.3719 

CP    * 

- 
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= 
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= 
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= 
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CP 
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= 

0.6750 

F3 
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= 
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F3 
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CP 

4c 

X 

0.6  800 

F3 

= 
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= 
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F  3 

- 
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CP 

* 

= 
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P3 

= 
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CP    * 

= 

0  .  6  90  0 

*  3 

- 

20  .7U6U 

CP 

4c 

= 
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F3 

= 
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CP    * 

= 
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F3 

- 
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CP 

4c 

= 
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F3 

= 
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= 
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F3 

= 
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CP 

4c 

= 
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-= 

21.6966 
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